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Atomic polarization phenomena impinge upon a number of areas and processes in physics. The 
dielectric constant and refractive index of any gas are examples of macroscopic properties that are 
largely determined by the dipole polarizability. When it comes to microscopic phenomena, the 
existence of alkaline-earth anions and the recently discovered ability of positrons to bind to many 
atoms are predominantly due to the polarization interaction. An imperfect knowledge of atomic 
polarizabilities is presently looming as the largest source of uncertainty in the new generation of 
optical frequency standards. Accurate polarizabilities for the group I and II atoms and ions of 
the periodic table have recently become available by a variety of techniques. These include refined 
many-body perturbation theory and coupled-cluster calculations sometimes combined with precise 
experimental data for selected transitions, microwave spectroscopy of Rydberg atoms and ions, 
refractive index measurements in microwave cavities, ab initio calculations of atomic structures 
using explicitly correlated wave functions, interferometry with atom beams, and velocity changes of 
laser cooled atoms induced by an electric field. This review examines existing theoretical methods 
of determining atomic and ionic polarizabilities, and discusses their relevance to various applications 
with particular emphasis on cold-atom physics and the metrology of atomic frequency standards. 

PACS numbers: 31.15.ap, 32.10.Dk, 42.50.Hz, 51.70.+f 



I. INTRODUCTION 

By the time Maxwell presented his article on a Dynam- 
ical Theory of the Electromagnetic Field it was un- 
derstood that bulk matter had a composition of particles 
of opposite electrical charge, and that an applied electric 
field would rearrange the distribution of those charges 
in an ordinary object. This rearrangement could be de- 
scribed accurately even without a detailed microscopic 
understanding of matter. For example, if a perfectly con- 
ducting sphere of radius ro is placed in a uniform electric 
field F, simple potential theory shows that the resulting 
electric field at a position r outside the sphere must be 
F — V(F • rrg/r 3 ). This is equivalent to replacing the 
sphere with a point electric dipole, 

d = aF, (1) 

where a — r$ is the dipole polarizability of the 
sphere 1 . An arbitrary applied electric field can be de- 
composed into multipole fields of the form F*(r) = 
—F k V{r k C k (v)), where C*(r) is a spherical tensor Q. 
Each of these will induce a multipole moment of F k r^ k+1 



in the conducting sphere, corresponding to a multipole 
polarizability of a k = r^ k+1 . Treatment of the electrical 
polarizabilities of macroscopic bodies is a standard topic 
of textbooks on electromagnetic theory, and the only ma- 
terial properties that it requires are dielectric constants 
and conductivities. 

Quantum mechanics, on the other hand, offers a fun- 
damental description of matter, incorporating the effects 
of electric and magnetic fields on its elementary con- 
stituents, and thus enables polarizabilities to be calcu- 
lated from first principles. The standard framework for 
such calculations, perturbation theory, was first laid out 
by Schrodinger 3] in a paper that reported his calcula- 
tions of the Stark effect in atomic hydrogen. A system of 
particles with positions and electric charges qi exposed 
to a uniform electric field, (F = FF), is described by the 
Hamiltonian 

H = H - FF ■ d, (2) 

where Hq is the Hamiltonian in the absence of the field, 
and d is the dipole moment operator, 

d = ^2q i r l . (3) 



For notational convenience, we use the Gaussian system of elec- 
trical units, as discussed in subsection A below. In the Gaussian 
system, electric polarizability has the dimensions of volume. 



Treating the field strength, F = |F|, as a perturbation 
parameter, means that the energy and wave function can 
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be expanded as 
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The first-order energy E\ = if |$o) is an eigenfunc- 
tion of the parity operator. In this case, |4>i) satisfies the 
equation 



(H -E )\$ 1 ) = -F-d\$ ) 



(6) 



From the solution to Eq. ([6]) , we can find the expectation 
value 



(tf|d|tf) = F(($ |d|$i) + ($i|d|$o)) 

= <5F (7) 

where a is a matrix. The second-order energy is given by 
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-F ■ aF. 
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(8) 



Although Eq. ^ can be solved directly, and in some 
cases in closed form, it is often more practical to express 
the solution in terms of the eigenfunctions and eigenval- 
ues of Hq, so that Eq. (jHJ takes the form 
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This sum over all stationary states shows that calculation 
of atomic polarizabilities is a demanding special case of 
the calculation of atomic structure. The sum extends in 
principle over the continuous spectrum, which sometimes 
makes substantial contributions to the polarizability. 

Interest in the subject of polarizabilities of atomic 
states has recently been elevated by the appreciation 
that the accuracy of next-generation atomic time and 
frequency standards, based on optical transitions 0-13, 
is significantly limited by the displacement of atomic en- 
ergy levels due to universal ambient thermal fluctuations 
of the electromagnetic field: blackbody radiation (BBR) 
shifts [Iol - fl3j . This phenomenon brings the most promis- 
ing approach to a more accurate definition of the unit of 
time, the second, into contact with deep understanding 
of the thermodynamics of the electromagnetic radiation 
field. 

Description of the interplay between these two fun- 
damental phenomena is a major focus of this review, 
which in earlier times might have seemed a pedestrian 
discourse on atomic polarizabilities. The precise calcu- 
lation of atomic polarizabilities also has implications for 
quantum information processing and optical cooling and 
trapping schemes. Modern requirements for precision 
and accuracy have elicited renewed attention to methods 



of accurate first-principles calculations of atomic struc- 
ture, which recently have been increased in scope and 
precision by developments in methodology, algorithms, 
and raw computational power. It is expected that the 
future will lead to an increased reliance on theoretical 
treatments to describe the details of atomic polarization. 
Indeed, at the present time, many of the best estimates of 
atomic polarizabilities are derived from a composite anal- 
ysis which integrates experimental measurements with 
first principles calculations of atomic properties. 

There have been a number of reviews and tabulations 
of atomic and ionic polarizabilities [l^ - [25| . Some of these 
reviews, e.g. 

[3 El HI El have lar § el y focussed upon 
experimental developments while others pjl [H[ have 
given theory more attention. 

In the present review, the strengths and limitations 
of different theoretical techniques are discussed in de- 
tail given their expected importance in the future. Dis- 
cussion of experimental work is mainly confined to pre- 
senting a compilation of existing results and very brief 
overviews of the various methods. The exception to this 
is the interpretation of resonance excitation Stark ion- 
ization spectroscopy [23[ since issues pertaining to the 
convergence of the perturbation analysis of the polariza- 
tion interaction are important here. The present review 
is confined to discussing the polarizabilities of low lying 
atomic and ionic states despite the existence of a body 
of research on Rydberg states [26|. High-order polariz- 
abilities are not considered except in those circumstances 
where they are specifically relevant to ordinary polariza- 
tion phenomenon. The influence of external electric fields 
on energy levels comprises part of this review as does the 
nature of the polarization interaction between charged 
particles with atoms and ions. The focus of this review 
is on developments related to contemporary topics such 
as development of optical frequency standards, quantum 
computing, and study of fundamental symmetries. Major 
emphasis of this review is to provide critically evaluated 
data on atomic polarizabilities. Table Q] summarizes the 
data presented in this review to facilitate the search for 
particular information. 



A. Systems of units 

Dipole polarizabilities are given in a variety of units, 
depending on the context in which they are determined. 
The most widely used unit for theoretical atomic physics 
is atomic units (a.u.), in which, e, m e , 4-7reo and the re- 
duced Planck constant h have the numerical value 1. 
The polarizability in a.u. has the dimension of vol- 
ume, and its numerical values presented here are thus 
expressed in units of ag, where ao ~ 0.052918 nm is the 
Bohr radius. The preferred unit systems for polarizabil- 
ities determined by experiment are A 3 , kHz/(kV/cm) 2 , 
cm 3 /mol or C-m 2 /N where C-m 2 /N is the SI unit. In 
this review, almost all polarizabilities are given in a.u. 
with uncertainties in the last digits (if appropriate) given 



3 



TABLE I: List of data tables. 



Table 


System 


Atoms and Ions 


States 


Data 


Table HV] 


Noble gases 


He, Ne, Ar, Kr, Xe, Rn, Li+, Na+, K+, Rb+, Cs+, Fr+ 
Be 2 +, Mg 2 +, Ca 2 +, Sr 2 +, Ba 2 +, Ra 2 + 


ground 


«o 


Table |V] 


Alkali atoms 


Li, Na, K, Rb, Cs, Fr 


ground ns, np 




Table EH 


Alkali ions 


Be+, Mg+, Ca+, Sr+, Ba+, Ra+ 


ground 


ao 


Table | VIII| 


Monovalent 


Li, Na, K, Ca+, Rb, Sr+ 


excited 


a , «2 


Table [VTTl 
Table E] 


Alkali atoms 
Alkali atom 


Resonance transition: Li, Na, K, Rb, Cs 
Na 


ns, np 
ground 


Aa 
ao 


Table |U 


Alkali atom 


Cs 


26 states 


ao, a2 


Table ED 


Group II type 


Be, Mg, Ca, Sr, Ba, Ra, A1+, Si 2+ , Zn, Cd, Hg, Yb 


ground, nsnp 3 Po 


ao 


Table IXHl 


Miscellaneous 


Al, Ga, In, Tl, Si, Sn, Pb, Ir, U, Cu, Ag, Au, 
A1+, Si 3+ , P 3+ , Kr 6+ , Cu+, Ag+, Hg+, Yb+, Zn+ 


ground 


ao 


Table |XIII| 


Miscellaneous 


Ca, Sr, Ba, Zn, Cd, Hg, Yb, Al, Tl, Yb+ 


excited 


ao, a2 


Table |XIV| 
Table |XVI| 

Table IXVIIII 
Table IXIXI 


Miscellaneous 
Miscellaneous 

Monovalent 
Alkali atoms 


Li, Na, Cs, Mg, Ca, Ba, Hg, Ga, Tl, Yb+ 
Mg, Ca, Sr, Yb, Zn, Cd, Hg, 
Ca+, Sr+, Hg+, Yb+, A1+, In+ 
Li, Na, K, Rb, Cs, Ba+, Yb+, Hg+ 
Li, Na, K, Rb, Cs, Fr 


clock 
transition 
ground hyperfine 
ground 


Aao 
Aj/bbr 

BBR 

c 6 



in parentheses. Conversion factors between the different 
units are listed in Table [TTJ The last line of the table 
gives conversion factors from SI units to the other units. 
For example, the atomic units for a can be converted to 
SI units by multiplying by 0.248832. 

Stark shift experiments which measure the change in 
photon frequency of an atomic transition as a function of 
electric field strength are usually reported as a Stark shift 
coefficient in units of kHz/(kV/cm) 2 . The polarizability 
difference is twice the size of the Stark shift coefficient, 
as in equation ([5]). 



II. ATOMIC POLARIZABILITIES AND 
FIELD-ATOM INTERACTIONS 



A. Static electric polarizabilities 

1. Definitions of scalar and tensor polarizabilities 

The overall change in energy of the atom can be evalu- 
ated within the framework of second-order perturbation 
theory. Upon reduction, the perturbation theory expres- 
sion given by Eq.® leads to a sum-over-states formula 
for the static scalar electric-dipole polarizability which is 
expressed most compactly in terms of oscillator strengths 



as 



a 



E 



fgn 



(AE n 



(10) 



In this expression, f gn is the absorption oscillator 
strength for a dipole transition from level g to level n, 
defined in a J-representation as 27] 



fgn — 



2|(V» g || rC'jr) || <P n )\ 2 AE ng 
3(2J„ + 1) 



(11) 



where A£L„ = E n — E g and C 1 (f ) is the spherical tensor 
of rank 1 [2j. The definition of the oscillator strength 
in LS coupling is transparently obtained from Eq. (fTTj) 
by replacing the total angular momentum by the orbital 
angular momentum. 

The polarizability for a state with non-zero angular 
momentum J depends on the magnetic projection M: 



3M 2 



a = ao + a 2 - 



J(J+1) 



J(2J- 1) 



(12) 



The quantity ao is called the scalar polarizability while 
«2 is the tensor polarizability in J representation. 

The scalar part of the polarizability can be determined 
using Eq. (|10[) . In terms of the reduced matrix elements 
of the electric-dipole operator, the scalar polarizability 
ao of an atom in a state ip with total angular momentum 
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TABLE II: Factors for converting polarizabilities between different unit systems. The table entries give the multiplying factor 
needed to convert the row entry to the corresponding column entry. The last column in the Table is the polarizability per mole 
and is often called the molar polarizability. The conversion factors from SI units to other units are given in the last line. Here, 
h is the Planck's constant, eo is the electric constant, ao is the Bohr radius, and JV4 is the Avogadro constant. 



A 



kHz/(kV/cm) 



C-m 2 /V 



cm 3 /mol 



A 3 

kHz/(kV/cm) 2 

Cm 2 /V 

cm 3 /mol 

Conversion from SI 



6.748335 
4.018778 
6.065100 xlO 4 
2.675205 
l/(47re ag) 



0.1481847 
1 

0.5955214 
..987552 xlO 39 
0.3964244 
10 30 /(47r eo ) 



0.2488319 
1.679201 
1 

6.626069 xl0~ 
0.6656762 
W~ 7 h 



1.648773 xl0~ 4i 
1.112650 xlO" 40 
1.509190 xlO 40 
1 

4.410816 xlO" 41 
1 



0.3738032 
2.522549 
1.502232 
2.267154 xlO 40 
1 

10 6 7v A /(3eo) 



J and energy E is also written as 

2 I^HrC 1 ^)!!^)! 2 



a 



3(2 J + 1) ^ E n — E 
The tensor polarizability a 2 is defined as 
5J(2J- 1) 



(13) 



a 2 = 4 



6(J+l)(2J+l)(2J + 3) 



1/2 



(14) 
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1 j 2 



E n — E 



It is useful in some cases to calculate polarizabilities in 
strict LS coupling. In such cases (28[, the tensor polar- 
izability «2,l f° r a state with orbital angular momentum 
L is given by 



a 2 , L = J2 




3(2£ + l) 



2|<y> |i rc 1 ^) || y> ra )|- 

E — E n 



(15) 



The tensor polarizabilities a 2 and a 2 ,L in the J and L 
representations, respectively, are related by 



a 2 



«2,i(-l) 



S+L+J+2 




For L = 1 and J = 3/2, Eq. (HU) gives a 2 = a 2iL . For 
L = 1. S = 1 and J = 1. Eq. (JIB) gives a 2 = -a 2yL /2. 
For i = 2, cv 2 = (7/10)a2,L for J = 3/2 and a 2 = a 2 ,L 
for J = 5/2. We use the shorter (?/>||jD||V> n ) designation 
for the reduced electric-dipole matrix elements instead of 
(V-llrCH?)!!^) below. 

Equation (|14l) indicates that spherically symmetric lev- 
els (such as the 6si/ 2 and 6pi/ 2 levels of cesium) only have 
a scalar polarizability. However, the hyperhne states of 



these levels can have polarizabilities that depend upon 
the hyperhne quantum numbers F and Mp. The rela- 
tionship between F and J polarizabilities is discussed in 
Ref. [29j. This issue is discussed in more detail in the 
section on BBR shifts. 

There are two distinctly different broad approaches to 
the calculation of atomic polarizabilities. The "sum-over- 
states" approach uses a straightforward interpretation of 
Eq. ((9l) with the contribution from each state ^ n being 
determined individually, either from a first principles cal- 
culation or from interpretation of experimental data. A 
second class of approaches solves inhomogeneous equa- 
tion ^ directly. We refer to this class of approaches as 
direct methods, but note that there are many different 
implementations of this strategy. 



2. The sum-over-states method 

The sum-over-states method utilizes expression such as 
Eqs. (JTU1 [131 - IT5]) to determine the polarizability. This 
approach is widely used for systems with one or two va- 
lence electrons since the polarizability is often dominated 
by transitions to a few low lying excited states. The sum- 
over-states approach can be used with oscillator strengths 
(or electric-dipole matrix elements) derived from exper- 
iment or atomic structure calculations. It is also possi- 
ble to insert high-precision experimental values of these 
quantities into an otherwise theoretical determination of 
the total polarizability. For such monovalent or diva- 
lent systems, it is computationally feasible to explicitly 
construct a set of intermediate states that is effectively 
complete. Such an approach is computationally more 
difficult to apply for atoms near the right hand side of 
the periodic table since the larger dimensions involved 
would preclude an explicit computation of the entire set 
of intermediate state wave functions. 

For monovalent atoms, it is convenient to separate the 
total polarizability of an atom into the core polarizabil- 
ity a corc and the valence part defined by Eq. (JT3) . The 
core contribution actually has two components, the po- 
larizability of the ionic core and a small change due to 
the presence of the valence electron [30]. For the alkali 
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atoms, the valence part of the ground state polarizabil- 
ity is completely dominated by the contribution from the 
lowest excited state. For example, the 5s — 5pi/2 and 
5s— 5p 3 / 2 transitions contribute more than 99% of the Rb 
valence polarizability [3l|. The Rb + core polarizability 
contributes about 3%. Therefore, precision experimental 
measurements of the transition rates for the dominant 
transitions can also be used to deduce accurate values 
of the ground state polarizability. However, this is not 
the case for some excited states where several transitions 
may have large contributions and continuum contribu- 
tions may be not negligible. 

This issue is illustrated using the polarizability of the 
5pi/2 state of the Rb atom [30j, which is given by 



ao(5pi/ 2 ) 



\(ns\\D\\5 Pl/2 )\ 2 

Ens ~~ E^Pi/2 



3 ^ 



1 V l(^3/2p||5pi/ 2 )| ; 



(17) 



3 4^ E nd 



'3/2 



E 



'5P1/2 



We present a solution to the Eq. p7|) that combines 
first principles calculations with experimental data. The 
strategy to produce a high-quality recommended value 
with this approach is to calculate as many terms as real- 
istic or feasible using the high-precision atomic structure 
methods. Where experimental high-precision data are 
available (for example, for the 5s — 5p transitions) they 
are used in place of theory, assuming that the expected 
theory uncertainty is higher than that of the experimen- 
tal values. The remainder that contains contributions 
from highly-excited states is generally evaluated using 
(Dirac-Hartree-Fock) DHF or random-phase approxima- 
tion (RPA) methods. In our example, the contribution 
from the very high discrete (n > 10) and continuum 
states is about 1.5% and cannot be omitted in a pre- 
cision calculation. Table IIIII lists the dipole matrix el- 
ements and energy differences required for evaluation of 
Eq. (fl7|) as well as the individual contributions to the 



polarizability. Experimental values from 32[ are used 
for the 5s — bpj matrix elements, otherwise the matrix 
elements are obtained from the all-order calculations of 
Ref. [13] described in Section IIVFI Absolute values of 
the matrix elements are given. Experimental energies 
from [33l l34j are used. Several transitions give signifi- 
cant contributions. This theoretical number agrees with 
experimental measurement within the uncertainty. The 
comparison with experiment is discussed in Section [Vj 



3. Direct methods 

From a conceptual viewpoint, the finite-field method 
represents one of the simplest ways to compute the po- 
larizability. In this approach, one solves the Schrodinger 
equation using standard techniques for the perturbed 
Hamiltonian given by Eq. ^ for a variety of values of F. 
The polarizability is then extracted from the dipole mo- 



TABLE III: The contributions (in a.u.) to the scalar polar- 
izability of the Rb atom in the 5pi/ 2 state [3(j. The uncer- 
tainties in each term are enclosed in parenthesis. The corre- 
sponding energy differences AE — E n — E$ Pl/2 [34| are given 
in cm -1 , which can be converted to atomic units by division 
by 219474.6. Experimental values from [321 ] are used for ab- 
solute values of the 5s — 5pj matrix elements, otherwise the 
matrix elements are obtained from all-order calculations of 
Ref. B. 



Contribution 


|<n||D||5p 1/2 )| 


AE 


Qo(5pi/ 2 ) 




- 5s 


4.231 


-12579 


-104.11(15) 


5pi/2 


- 6s 


4.146 


7554 


166.5(2.2) 


5pi/2 


- 7s 


0.953 


13733 


4.835(16) 


5pi/2 


-8s 


0.502 


16468 


1.120(7) 


5P1/2 


-9s 


0.331 


17920 


0.448(3) 


5P1/2 


- 10s 




±0 ( 00 




5P1/2 


- lis 


0.189 


19338 


0.135(1) 


5P1/2 


- (12 - oo)s 






1.9(0.2) 


5P1/2 


- 4d 3/2 


8.017 


6777 


694(30) 


5P1/2 


- 5d 3 / 2 


1.352 


13122 


10.2(9) 


5P1/2 


- 6d 3 / 2 


1.067 


16108 


5.2(1.1) 


5P1/2 


- 7d 3/2 


0.787 


17701 


2.6(4) 


5P1/2 


— 8d 3 / 2 


0.605 


18643 


1.4(2) 


5P1/2 


- 9d 3/2 


0.483 


19243 


0.89(10) 


5P1/2 


- (10 - oo)d 3/2 






10.5(10.5) 










9.08(45) 


Total 








805(31) 



ment or the energy eigenvalues of the perturbed Hamil- 
tonian. This usually entails doing a number of calcula- 
tions at different discrete field strengths. This approach 
is generally used to obtain polarizabilities in coupled- 
cluster calculations (see, for example, Refs. [35, 36]). We 
note that linearized coupled-cluster calculations are im- 
plemented very differently, and sum-over-states is used 
for the polarizability calculations [37J. These differences 
between coupled-cluster calculations are discussed in Sec- 
tion [Tv] 

Another direct approach to calculating polarizabil- 
ity is the perturbation- variation method [381 ] . The 
perturbation- variation approach has been outlined in the 
introduction as Eqs. §5§ to (J7J. The unperturbed state, 
|$o) an d perturbed state, |<I>i) would be written as a 
linear combinations of basis states. Equations ([6]) and 
(|7|) then reduce to sets of matrix equations. A general 
technique for solving the inhomogeneous equation ^ has 
been described by Dalgarno and Lewis in Ref. [39j | . 

Exact solutions to Eqs. ([5]) - ([7]) are possible for atomic 
hydrogen and hydrogenic ions. The non-relativistic so- 
lutions were first obtained independently in 1926 by Ep- 
stein [13], Waller [4l|, and Wentzel [J|; the relativis- 
tic case remains a subject of current research interest 
[43M46| . The nonrelativistic equations are separable in 
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parabolic coordinates, and the polarizability of a hydro- 
genic ion of nuclear charge Z in the state \nin 2 m > is (in 
a.u.) 



3(ni - n 2 ) 2 - 9m 2 + 19]og, (18) 



where n\ , n 2 are parabolic quantum numbers [47] , m is 
the projection of the orbital angular momentum onto the 
direction of the electric field, and n = tlx + n 2 + |m| + 1 
is the principal quantum number. A convenient special 
case is n = \m\ + 1, which corresponds to the familiar 
circular states of hydrogen in spherical coordinates, with 
orbital angular momentum I = \m\ = n — 1; for these 
states, a = (|m| + 2)(|m| + 9/4)ag. 

For the H Is ground state exposed to an electric field 
F = Fz, the solution to Eq. © is (in a.u.) 



(19) 
(20) 



$o = e 7Vtt, 

$i = -*(l + r/2)$ , 



from which a — (9/2)oq. Note that although $i of 
Eq. (|20p is a p state, it is much more compact than 
any of the discrete np eigenstates of H. Thus building 
up $i by the sum-over-states approach requires a signif- 
icant contribution from the continuous spectrum of H. 
This is depicted in Fig. [TJ which employs the histogram 
construction of Fano and Cooper [48| to show the con- 
nection between discrete and continuum contributions to 
the sum over states. About 20% of the polarizability of 
H Is comes from the continuum. 



FIG. 1: Solid line: histogram representation of the sum-over- 
states contributions to the polarizability of H Is. Following 
Ref. [13], the contributions of discrete states (e.g. 2p) are 
spread over the inverse density of states, to show continuity 
with the continuum contributions near energy E = 0. The 
polarizability, a, is equal to the area under this curve. Dashed 
line: the same construction, for an electron bound to a one- 
dimensional delta- function potential with energy E = —1/2 
a.u. From [ig|| . 

Clearly, the direct solution of the Schrodinger equa- 
tion for an atom in the presence of an electric field and 



subsequent determination of the polarizability is formally 
equivalent to the sum-over-states approach described in 
the previous subsection. However, it is useful to com- 
ment on how this equivalence is actually seen in calcu- 
lations for many-electron atoms. For example, random- 
phase-approximation (RPA) results for polarizabilities of 
closed-shell atoms [50] that were obtained by direct so- 
lution of inhomogeneous equation are the same (up to 
numerical uncertainty of the calculations) as sum-over- 
state RPA results obtained using formula 



\m\D^u m )^ a \\D^u m )\ 



3 ^ 

ma 



E,, 



E a 



(21) 



where (ip a \\D \\ip m ) is reduced matrix element of 
dipole operator obtained in the DHF approximation and 
the (?/> a ||-D RPA ||?/' m ) matrix elements include RPA terms 
using many-body perturbation theory as discussed, for 
example, in [HTJ] . The index a refers to all core orbitals, 
while the m includes all other orbitals. The sum-over- 
states can be calculated with a finite basis set [Hj], and 
such an approach intrinsically includes the continuum 
states when complete sum over the entire basis set is 
carried out. When the contributions from highly-excited 
states are significant, it becomes difficult to account for 
these terms accurately within the framework of the sum- 
over-states approach. Direct method automatically ac- 
counts for these states and this problem does not arise. 
However, it becomes difficult and cumbersome to include 
corrections to the dipole operator beyond RPA. The 
method implemented in [Hoi ] is different from the finite 
field approach and does not involve performing a num- 
ber of calculations at different discrete field strengths. 

In most high-precision calculations, the determina- 
tion of polarizabilities follows the calculation of wave 
functions or quantities that represent the wave func- 
tions (such as excitation coefficients). The type of ap- 
proach used for this initial calculation generally deter- 
mines whether polarizabilities are determined by Eqs. ([5]) 
or by sum-over-states method. For example, relativistic 
linearized coupled-cluster approach (37j is formulated in 
a way that does not explicitly generate numerical wave 
functions on a radial grid, and all quantities are expressed 
in terms of excitation coefficients. Therefore, the polar- 
izabilities are calculated by the sum-over-states method 
using resulting high-quality dipole matrix elements and 
energies. In the case of methods that combine rela- 
tivistic configuration interaction and perturbation the- 
ory [CI+MBPT], it is natural to determine polarizabili- 
ties by directly solving the inhomogeneous equation. In 
this case, it is solved in the valence space with the ef- 
fective operators that are determined using MBPT [53j . 
The ionic core polarizability is calculated separately in 
this approach. The effective dipole operator generally 
includes RPA corrections, with other corrections calcu- 
lated independently. 

The direct and sum-over-states approaches can also be 
merged in a hybrid approach. One strategy is to perform 
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a direct calculation using the best available techniques, 
and then replace the transition matrix elements for the 
most important low-lying states with those from a higher 
level theory. This hybrid method is discussed further in 
the sections on the CI+MBPT and CI+all-order meth- 
ods. 



B. The frequency-dependent polarizability 

So far, we have described the polarizability for static 
fields. The numerical value of the polarizability changes 
when the atom is immersed in an alternating (AC) elec- 
tromagnetic field. To second order, one writes AE = 
— ^a(oj)F 2 + . . .. The valence part of the scalar fre- 
quency dependent polarizability, usually called the dy- 
namic polarizability, is calculated using the sum-over- 
states approach with a straightforwardly modified ver- 
sion of Eq. (TTT 



a (ui) 



3(2 J + 1) 



AEM\DW n )f 
(AE) 2 - w 2 



(22) 



Eq. (|22l) assumes that oj is at least a few linewidths away 
from resonant frequencies defined by AE = E n — E. As 
noted previously, atomic units are used throughout this 
paper, and h = 1. The core part of the polarizability may 
also be corrected for frequency dependence in random 
phase approximation by similarly modifying the formula 
(f2"Tj) . Static values may be used for the core contribu- 
tion in many applications since the frequencies of inter- 
est, (i.e. corresponding to commonly used lasers) are very 
far from the excitation energies of the core states. The 
calculations of the ground and excited state frequency- 
dependent polarizabilities of the alkali-metal atoms are 
described in detail in Refs. [54| and [29j, respectively. It 
is essentially the same as the calculation of the static 
polarizability described in Section Hi A 21 only for uj ^ 0. 

The expression for the tensor polarizability given by 
Eq. (114[) is modified in the same way, i.e. by replacing 



1 AE 

AE ^ AE 2 — i 



(23) 



There has been more interest recently in the de- 
termination of frequency-dependent polarizabilities due 
to the need to know various "magic wavelengths" [55[ 
for the development of optical frequency standards and 
other applications. At such wavelengths, the frequency- 
dependent polarizabilities of two states are the same, and 
the AC Stark shift of the transition frequency between 
these two states is zero. An example of the calculation 
of frequency-dependent polarizabilities and magic wave- 
lengths is given in Section [VII Bl Experimentally deter- 
mined magic wavelengths may also be used to gauge the 
accuracy of the theory. 



III. MEASUREMENTS OF POLARIZABILITIES 
AND RELATED QUANTITIES 

Experimental measurements of atomic and ionic polar- 
izabilities are somewhat rarer than theoretical determina- 
tions. There are two types of measurements, those which 
directly determine the polarizability, and those which de- 
termine differences in polarizabilities of two states from 
Stark shift of atomic transitions. 

For the most part, we make brief comments on the 
major experimental techniques and refer the reader to 
primarily experimental reviews [l7l [l9l [22l [23j j for further 
details. 



A. /-sum rules 



This approach makes use of Eqs. (| 10H1 5|) . Many of the 
most interesting atoms used in cold atom physics typi- 
cally have only one or two valence electrons. The ground 
state polarizability of these atoms is dominated by a sin- 
gle low- lying transition. As mentioned in Section [II A 21 
97% of the total value of Rb ground state polarizability 
comes from 5s — > bp transition. In the case of Na, about 
99.4% of the valence polarizability and 98.8% of the total 
polarizability of sodium arises from the 3s — » 3p resonant 
transition. 

Composite estimates of the polarizability using both 
experimental and theoretical inputs are possible. One 
type of estimate would use experimental oscillator 
strengths to determine the valence polarizability. This 
could be combined with a core contribution obtained by 
other methods to estimate the total polarizability. An- 
other approach replaces the most important matrix ele- 
ments in a first-principles calculation by high precision 
experimental values [5a . [52} ■ Various types of experi- 
ments may be used to determine particular matrix ele- 
ments, including photo-association experiments [58J, life- 
time, oscillator strengths, or Stark shift measurements 
[30| with photoassociation experiments generally giving 
the most reliable matrix elements. This hybrid method 
may provide values accurate to better than 0.5% in cer- 
tain cases 15611 . 



B. Dielectric constant 

The dielectric constant K of an atomic or molecular 
gas is related to the dipole polarizability, a, by the iden- 
tity 



a 



K - 1 
AttN ' 



(24) 



where N is the atomic number density. The technique 
has only been applied to the rare gas atoms, and the 
nitrogen and oxygen atoms by the use of a shock tube. 
Results for the rare gases typically achieve precisions of 
0.01 — 0.1%. Examples are reported in Table ITVl 
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C. Refractive index 

The frequency-dependent refractive index of a gas 
n(w), is related to the polarizability by the expression 



a(u>) 



n(uj) — 1 
2ttN ' 



(25) 



where N is the atomic number density. The static dipole 
polarizability, a(0), can be extracted from the frequency- 
dependent polarizability a{uj) by the following technique. 

The energy denominator in Eq. (|22|) can be expanded 
when the frequency is smaller than the frequency of the 
first excitation giving 



o(w) = a(0) + w 2 5(-4) + w 4 5(-6) . 



(26) 



The S(—q) factors are the Cauchy moments of the oscil- 
lator strength distribution and are defined by 



fgn 



(27) 



Specific Cauchy moments arise in a number of atomic 
physics applications, as reviewed by Fano and Cooper 
[48|. For example, the Thomas- Reiche-Kuhn sum rule 
states that 5(0) is equal to the number of electrons in the 
atom. The 5(— 3) moment is related to the non-adiabatic 
dipole polarizability [53, [6(| ■ 

The general functional dependence of the polarizability 
at low frequencies is given by Eq. (|2"6"j) (6ll l62l |. The 
achievable precision for the rare gases is 0.1% or better 
[HI, [63| . Experiments on the vapours of Zn, Cd and Hg 
gave polarizabilities with uncertainties of 1 — 10% [63 . [65| . 



D. Deflection of an atom beam by electric fields 

The beam deflection experiment is conceptually sim- 
ple. A collimated atomic beam is directed through an 
interaction region containing an inhomogeneous electric 
field. While the atom is in the interaction region, the 
electric field F induces a dipole moment on the atom. 
Since the field is not uniform, a force proportional to the 
gradient of the electric field and the induced dipole mo- 
ment results in the deflection of the atomic beam. The 
polarizability is deduced from the deflection of the beam. 
The overall uncertainty in the derived polarizabilities is 
between 5 — 10% (66|. Therefore, this method is mainly 
useful at the present time for polarizability measurements 
in atoms inaccessible by any other means. 



giving a magnetic deflection force in addition to the elec- 
tric deflection. The experiment is tuned so that the elec- 
tric and magnetic forces are in balance. The polarizabil- 
ity can be determined since the magnetic moments of 
many atoms are known. Uncertainties range from 2% to 
10% MM- 



F. Atom interferometry 

The interferometry approach splits the beam of atoms 
so that one path sends a beam through a parallel plate 
capacitor while the other goes through a field free re- 
gion. An interference pattern is then measured when the 
beams are subsequently merged and detected. The po- 
larizability is deduced from the phase shift of the beam 
passing through the field free region. So far, this ap- 
proach has been used to measure the polarizabilities of 
helium (see [70| for a discussion of this measurement), 
lithium [7l[ , sodium [72|, [73[ , potassium [73| , and rubid- 
ium ill achieving uncertainties of 0.35 — 0.8%. 

It has been suggested that multi-species interferom- 
eters could possibly determine the polarizability ratio 
R = ax /ay to 10~ 4 relative accuracy |70j |. Conse- 
quently, a measurement of R in conjunction with a known 
standard, say lithium, could lead to a new level of pre- 
cision in polarizability measurements. Already the Na:K 
and Na:Rb polarizability ratios have been measured with 
a precision of 0.3% [73j |. 



G. Cold atom velocity change 

The experiment of Amini and Gould [74j measured the 
kinetic energy gained as cold cesium atoms were launched 
from a magneto-optical trap into a region with a finite 
electric field. The kinetic energy gained only depends 
on the final value of the electric field. The experimen- 
tal arrangement actually measures the time of return for 
cesium atoms to fall back after they are launched into a 
region between a set of parallel electric-field plates. The 
only such experiment reported so far gave a very precise 
estimate of the Cs ground state polarizability, namely 
ao — 401.0(6) a.u.. This approach can in principle be ap- 
plied to measure the polarizability of many other atoms 
with a precision approaching 0.1% (22j. 



H. Other approaches 



E. The E-H balance method 

In this approach, the E-H balance configuration ap- 
plies an inhomogeneous electric field and an inhomoge- 
neous magnetic field in the interaction region [67). The 
magnetic field acts on the magnetic moment of the atom 



The deflection of an atomic beam by pulsed lasers has 
been used to obtain the dynamic polarizabilities of ru- 
bidium and uranium (75l . l76j . The dynamic polarizabili- 
ties of some metal atoms sourced from an exploding wire 
have been measured interferometrically [ill . These 
approaches measure polarizabilities to an accuracy of 5- 
20%. 
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I. Spectral analysis for ion polarizabilities 

The polarizability of an ion can in principle be ex- 
tracted from the energies of non-penetrating Rydberg se- 
ries of the parent systemlU [j^, [8(| ■ The polarizability 
of the ionic core leads to a shift in the (n, L) energy levels 
away from their hydrogenic values. 

Consider a charged particle interacting with an atom 
or ion at large distances. To zeroth order, the interac- 
tion potential between a highly excited electron and the 
residual ion is just 



V(r) = 



Z - N 



(28) 



where Z is the nuclear charge and N is the number 
of electrons. However, the outer electron perturbs the 
atomic charge distribution. This polarization of the 
electron charge cloud leads to an attractive polarization 
potential between the external electron and the atom. 
The Coulomb interaction in a multipole expansion with 
I r |>| x I, is written as 



The non-adiabatic interaction is repulsive for atoms in 
their ground states. The polarization interaction includes 
further adiabatic, non-adiabatic and higher order terms 
that contribute at the r~ 7 and r -8 , but there has been 
no systematic study of what could be referred to as the 
non-adiabatic expansion of the polarization potential. 

When the Rydberg electron is in a state that has neg- 
ligible overlap with the core (this is best achieved with 
the electron in high angular momentum orbitals), then 
the polarization interaction usually provides the domi- 
nant contribution to this energy shift. Suppose the dom- 
inant perturbation to the long-range atomic interaction 
is 



- C4 v i , 

V P o\{r) = -— - — 



r 6 ■ 



(34) 



where C 4 = a /2 and C 6 = (a - 6/3)/2. Equation ([55]) 
omits the CV / r 7 and C% /r 8 terms that are included in a 
more complete description (82H84j . The energy shift due 
to an interaction of this type can be written 



1 



= E Cfe (x)-C*(r)-^ T . 



(29) 



AE 
A(F^ 



= C4 + C{ 



--6\ 



A(r- 4 ) 



(35) 



Applying second-order perturbation theory leads to the 
adiabatic polarization potential between the charged par- 
ticle and the atom, e.g. 



^poi(r) 



00 



,,-Efe 



k=l 



(30) 



The quantities a Ek are the multipole polarizabilities de- 
fined as 



Rk 



E 



/.in 

Jgn 



(k) 



{AE gn f 



(31) 



In this notation, the electric-dipole polarizability is 
written as a E1 , and fg„ is the absorption oscillator 
strength for a multipole transition from g — > n. Equa- 
tion (|3"0"|) . with its leading term involving the dipole po- 
larizability is not absolutely convergent in k [81j. At 
any finite r, continued summation of the series given by 
Eq. (|30|) . with respect to k, will eventually result in a 
divergence in the value of the polarization potential. 

Equation (|30[) is modified by non-adiabatic corrections 
59, 601. The non-adiabatic dipole term is written as 



non — ad 



6/3o 
2r 6 



(32) 



where the non-adiabatic dipole polarizability, (3q is de- 
fined 



f (i) 

Jgn 



^ 2(A£ sn )3 



(33) 



where AE is usually the energy difference between two 
Rydberg states. The expectation values A(r~ 6 ) and 
A(r~ 4 ) are simply the differences in the radial expecta- 
tions of the two states. These are easily evaluated using 
the identities of Bockasten [85| . Plotting 



A(r-*> 



AE 

A(r-*> 



versus 



yields C4 as the intercept and Cq as the gradient. 
Such a graph is sometimes called a polarization plot. 

Traditional spectroscopies such as discharges or laser 
excitation find it difficult to excite atoms into Rydberg 
states with L > 6. Exciting atoms into states with L > 6 
is best done with resonant excitation Stark ionization 
spectroscopy (RESIS) [23(. RESIS spectroscopy first ex- 
cites an atomic or ionic beam into a highly-excited state, 
and then uses a laser to excite the system into a very 
highly-excited state which is Stark ionized. 

While this approach to extracting polarizabilities from 
Rydberg series energy shifts is appealing, there are a 
number of perturbations that act to complicate the 
analysis. These include relativistic effects AE rc \ : Stark 
shifts from ambient electric fields AE SS , second-order ef- 
fects due to relaxation of the Rydberg electron in the 
field of the polarization potential AE SCC [86-^88j, and fi- 
nally the corrections due to the Cr/r 7 and Cs/r 8 terms, 
AE 7 ,AE$, and AE 8L . Therefore, the energy shift be- 
tween two neighbouring Rydberg states is 

AE = AE 4 + AE 6 + AE 7 + AE 8 + AE SL 

+ AE Icl + AE SCC + AE SS . (36) 

One way to solve the problem is to simply subtract these 
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terms from the observed energy shift, e.g. 

AE C = AE ohs 
A(r- 4 ) ~ A(?- 4 > 

{ AE rcl + AE SCC + AE S 



\ A(r- 4 ) 

/ AE7 + AE$ + AE$l 

{ A(r-4) 



(37) 



and then deduce C4 and Cq from the polarization plot of 
the corrected energy levels [84[ . 



Stark shift measurements of polarizability 
differences 



The Stark shift experiment predates the formulation of 
quantum mechanics in its modern form (89J. An atom is 
immersed in an electric field, and the shift in wavelength 
of one of its spectral lines is measured as a function of the 
field strength. Stark shift experiments effectively mea- 
sure the difference between the polarizability of the two 
atomic states involved in the transition. Stark shifts can 
be measured for both static and dynamic electric fields. 
While there have been many Stark shift measurements, 
relatively few have achieved an overall precision of 1% or 
better. 

While the polarizabilities can generally be extracted 
from the Stark shift measurement, it is useful to compare 
the experimental values directly with theoretical predic- 
tions where high precision is achieved for both theory and 
experiment. In this review, comparisons of the theoreti- 
cal static polarizability differences for the resonance tran- 
sitions involving the alkali atoms with the corresponding 
Stark shifts are provided in Section [V] Some of the alkali 
atom ex per iments report precisions between 0.01 and 0.1 
a.u. [90j-|93| . The many Stark shift experiments involving 
Rydberg atoms [94| are not detailed here. 

Selected Stark shifts for some non-alkali atoms that 
are of interest for applications described in this review 
are discussed in Section [V] as well. The list is restricted 
to low-lying excited states for which high precision Stark 
shifts are available. When compared with the alkali 
atoms, there are not that many measurements and those 
that have been performed have larger uncertainties. 

The tensor polarizability of an open shell atom can 
be extracted from the difference in polarizabilities be- 
tween the different magnetic sub-levels. Consequently, 
tensor polarizabilities do not rely on absolute polariz- 
ability measurements and can be extracted from Stark 
shift measurements by tuning the polarization of a probe 
laser. Tensor polarizabilities for a number of states of 
selected systems are discussed in Section IVl 

One unusual experiment was a measurement of the AC 
energy shift ratio for the 6s and 5d 3 / 2 states of Ba + to 



K. AC Stark shift measurements 

There are few experimental measurements of AC Stark 
shifts at optical frequencies. Two recent examples would 
be the determination of the Stark shift for the Al + clock 
transition [97j and the Li 2s-3s Stark shifts 98 1 at the 
frequencies of the pump and probe laser of a two-photon 
resonance transition between the two states. One diffi- 
culty in the interpretation of AC Stark shift experiments 
is the lack of precise knowledge about the overlap of the 
laser beam with atoms in the interaction region. This 
is also a complication in the analysis of experiment on 
deflection of atomic beams by lasers [75|, [76[ . 



IV. PRACTICAL CALCULATION OF ATOMIC 
POLARIZABILITIES 

There have been numerous theoretical studies of 
atomic and ionic polarizabilities in the last several 
decades. Most methods used to determine atomic wave 
functions and energy levels can be adapted to generate 
polarizabilities. These have been divided into a number 
of different classes that are listed below. We give a brief 
description of each approach. It should be noted that the 
list is not exhaustive, and the emphasis here has been on 
those methods that have achieved the highest accuracy 
or those methods that have been applied to a number of 
different atoms and ions. 



A. Configuration interaction 

The configuration interaction (CI) method [9!| and its 
variants are widely used for atomic structure calculations 
owing to general applicability of the CI method. The 
CI wave function is written as a linear combination of 
configuration state functions 



(38) 



an accuracy of 0.11% [95[. This experiment does not give 
polarizabilities, and is mainly valuable as an additional 
constraint upon calculation 96]. 



i.e. a linear co mbina tion of Slater determinants from a 
model subspace [l OOj j . Each configuration is constructed 
with consideration given to anti-symmetrization, angu- 
lar momentum and parity requirements. There is a great 
deal of variety in how the CI approach is implemented. 
For example, sometimes the exact functional form of the 
orbitals in the excitation space is generated iteratively 
during successive diagonalization of the excitation basis. 
Such a scheme is called the multi-configuration Hartree- 
Fock (MCHF) or mul ti-co nfiguration self consistent field 
(MCSCF) approach [lOlj . The relativistic version of 
MCHF is referre d to as multi-configuration Dirac-Fock 
(MCDF) method flO^ . 

The CI approach has a great deal of generality since 
there are no restrictions imposed upon the virtual or- 
bital space and classes of excitations beyond those limited 
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by the computer resources. The method is particularly 
useful for open shell systems which contain a number of 
strongly interacting configurations. On the other hand, 
there can be a good deal of variation in quality between 
different CI calculations for the same system, because 
of the flexibility of introducing additional configuration 
state functions. 

The most straightforward way to evaluate polarizabil- 
ity within the framework of the CI method it to use 
a direct approach by solving the inhomogeneous equa- 
tion ([5]). RPA corrections to the dipole operator can be 
incorporated using the effective operator technique de- 
scribed in Section IIVG1 It is also possible to use CI- 
generated matrix elements and energies to evaluate sums 
over states. The main drawback of the CI method is its 
loss of accuracy for heavier systems. It becomes difficult 
to include a sufficient number of configurations for heav- 
ier systems to produce accurate results even with modern 
computer facilities. One solution of this problem is to use 
a semi-empirical core potential (CICP method) described 
in the next subsection. Another, ab initio solution, in- 
volves construction of the effective Hamiltonian using ei- 
ther many-body perturbation theory (CI+MBPT) or all- 
order linearized coupled-cluster method (CI+all-order) 
and carrying out CI calculations in the valence sector. 
These approaches are described in the last two sections 
of this chapter. 



B. CI calculations with a semi-empirical core 
potential (CICP) 

The ab initio treatment of core-valence correlations 
greatly increases the complexity of any structure calcula- 
tion. Consequently, to include this physics in the calcu- 
lation, using a semi-empirical approach is an attr activ e 
alternative for an atom with a few valence electrons 

[IM- 

fl05j . 

In this method, the active Hamiltonian for a system 
with two valence electrons is written as 



H = ^2(-^V 2 i +V dh (r i ) + V exc (r i ) + V pl (r i ) 



— + V p2 (r 1 ,r 2 ). 
r\2 



(39) 



The Vdir and V^ xc represent the direct and exchange 
interactions with the core electrons. In some appro aches , 
thes e terms are represented by model potentials, |l06l - 
1 1Q8( | . More refined approaches evaluate Vdir and y exc us- 
ing core wave functions c alcul a ted with the Hartree-Fock 
(or Dirac-Fock) method [IM EM EH • The one-body 
polarization interaction V po i(r) is semi-empirical in na- 
ture and can be written in its most general form as an 
^-dependent potential, e.g. 



^Pi(r) 



where a is the static dipole polarizability of the core and 
gj (r) is a cutoff function that eliminates the 1/r 4 singu- 
larity at the origin. The cutoff functions usually include 
an adjustable parameter that is tuned to reproduce the 
binding energies of the valence states. The two-electron 
or di-electronic polarization potential is written 



Vp 2 (r i ,r :j ) 



/>»3 (yi 3 

i 3 



(41) 



There is variation between expressions for the core po- 
larization potential, but what is described above is fairly 
representative. One cho ice for the cutoff function is 
g\(r) = 1 — exp(— r 6 //^j) |l05j . but other choices exist. 

A complete treatment of the core-polarization correc- 
tions also implies that corrections have to be made to 
the multipole operators [IM EM, EH EH| ■ The modi- 
fied transition operator is obtained from the mapping 



£ C fc (r) ge (r)r k C k (r) 



(42) 



i m 



2r 4 



(40) 



Usage of the modified operator is essential to the correct 
prediction of the oscillator strengths. For exam ple, it 
reduces the K(4s Ap) oscillator strength by 8% [l04| . 

One advantage of this configuration interaction plus 
core-polarization (CICP) approach is in reducing the size 
of the calculation. The elimination of the core from ac- 
tive consideration permits very accurate solutions of the 
Schrodinger equation for the valence electrons. Intro- 
duction of the core-polarization potentials, V p \ and V p2 , 
introduces an additional source of uncertainty into the 
calculation. However, this additional small source of un- 
certainty is justified by the almost complete elimination 
of computational uncertainty in the solution of the re- 
sulting simplified Schrodinger equation. 

The CICP approach only gives the polarizability of the 
valence electrons. Core polarizabilities are typically quite 
small for the group I and II atoms, e.g. the ces ium atom 
has a large core polarizability of about 15.6 Oq [l 12l | . but 
this represents only 4% of the total ground state polariz- 
ability of 401 a,Q [74|. Hence, usage of moderate accuracy 
core polarizabilities sourced from theory or experiment 
will lead to only small inaccuracies in the total polariz- 
ability. 

Most implementation of the CICP approach to the 
calculation of polarizabilities have been within a non- 
relativistic framework. A relativistic variant (RCICP) 
has r ecently been applied to zinc, cadmium, and mercury 
[l 13| . It should be noted that even non- relativistic cal- 
culations incorporate relativistic effects to some extent. 
Tuning the core polarization correction to reproduce the 
experimental binding energy partially incorporates rela- 
tivistic effects on the wave function. 



C. Density functional theory 

Approaches based on Density Functional Theory 
(DFT) are not expected to give polarizabilities as accu- 



12 



rate as those coming from the refined ab initio calcula- 
tions described in the following sections. Polarizabilities 
from DFT calculations are most likely to be useful for 
systems for which large scale ab initio calculations are 
difficult, e.g. the transition metals. DFT calculations 
are often much less computationally expensive than ab 
initio calculations. There have bee n two relatively ex- 
tensive DFT compilations |ll4 , Ilia ] that have reported 
dipole polarizabilities for many atoms in the periodic ta- 
ble. 



D. Correlated basis functions 

The accuracy of atomic structure calculations can be 
dramatically improved by the use of basis functions which 
explicitly include the electron-electron coordinate. The 
most accurate calculations reported for atoms and ions 
with two or three electrons have typically been per- 
formed with exponential basis functions including the 
inter-electronic coordinates as a linear factor. A typical 
Hylleraas basis function for lithium would be 



X 



Jl ~J2 „33 JJX1 „J13 „J23 

'1 '2 '3 ' 



12 r{f r% exp (-an - /3r 2 - 7 r 3 ) . (43) 



Difficulties with performing the multi-center integrals 
have effectively precluded the use of such basis func- 
tions for systems with more than three electrons. Within 
the framework of the non-relativistic Schrodinger equa- 
tion, calculations with Hyllera as ba sis sets achieve accu- 
racies of 13 significant digits [l!6j | for polarizability of 
two-electron systems and 6 significa nt digits for the po- 
larizability of three-electron systems |ll7l . lll8| . Inclusion 
of relativistic and quantum electrodynamic (QED) cor- 
rections to t he po l ariza bility of helium has been carried 
out in Refs. [l!6l Ill9j |. and the resulting final value is 
accurate to 7 significant digits. 

Another correlated basis set that has recently found 
increasingly widespread use utilizes the explicitly corre- 
lated gaussian (ECG). A typical spherically symmetric 
explicitly cor relat ed gaussian for a three-electron system 
is written as Il2dll 



(44) 



The multi-center integrals that occur in the evaluation 
of the Hamiltonian can be generally reduced to analytic 
expressions that are relatively easy to compute. Cal- 
culations using correlated gaussians do not achieve the 
same precision as Hylleraas forms, but are still capable 
of achieving much higher precision than orbital based 
calculation s provide d the parameters ai and fiij are well 
optimized [lil [ml. 



E. Many-body perturbation theory 

The application of many-body perturbation theory 
(MBPT) is discussed in this section in the context of the 
Dirac equation. While MBPT has been applied with the 
non-relativistic Schrodinger equation, many recent ap- 
plications most relevant to this review have been using a 
relativistic Hamiltonian. 

The point of departure for the discussions of relativis- 
tic many-body perturbation theory (RMBPT) calcula- 
tions is the no-pair Hamiltonian obtained from QED by 
Brown and Ravenhall [l22j , where the contributions from 
negative-energy (positron) states are projected out. The 
no-pair Hamiltonian can be written in second-quantized 
form as H = Hq + V, where 



H = y^ei[alai] , 



V 



- (9ijkl + bijki) [a\a]aia k } 



(45) 
(46) 



ijkl 



Y (^dhf + Bdhf - U) l3 [a\a s ], 



and a c-number term that just provides an additive con- 
stant to the energy of the atom has been omitted. 

In Eqs. (03] -US]), a\ and are creation and annihila- 
tion operators for an electron state i, and the summation 
indices range over electron bound and scattering states 
only. Products of operators enclosed in brackets, such 
as [ajajajafe], designate normal products with respect to 
a closed core. The core DHF potential is designated by 
Vdhf and its Breit counterpart is designated by £?dhf- 
The quantity €i in Eq. (|45|) is the eigenvalue of the Dirac 
equation. The quantities gijki and bijki in Eq. (|46|) are 
two-electron Coulomb and Breit matrix elements, respec- 
tively 



9ijkl 

hjki 



1 



>'12 



MJ , (47) 
a.i-OL 2 + (e*i • r 12 )(£*2 ■ rn) 



2ru 



£7 



(48) 



where a. are Dirac matrices. 

For neutral atoms, the Breit interaction is often a 
small perturbation that can be ignored compared to the 
Coulomb interaction. In such cases, it is particularly con- 
venient to choose the starting potential U(r) to be the 
core DHF potential U — Vdhf, 



(^Hf)^ — Y^ \diaja ~ 9iaaj] 



(49) 



since with this choice, the second term in Eq. (|46|) van- 
ishes. The index a refers to all core orbitals. The Breit 
(£?DHF)r/ term is defined as 



(-BDHF)y — J]] [biaja ~ bi aa j] 



(50) 
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For monovalent atoms, the lowest-order wave function 
is written as 

l*i 0) )=4|0c), (51) 

where |0 C ) = a\al ■ ■ ■ aJjO) is the closed core wave func- 
tion, |0) being the vacuum wave function, and a \ being a 
valence-state creation operator. The indices a and b refer 
to core orbitals. 

The perturbation expansion for the wave function leads 
immediately to a perturbation expansion for matrix ele- 
ments. Thus, for the one-particle operator written in the 
second-quantized form as 

Z = ^2,Zija\a h (52) 

ij 

perturbation theory leads to an order-by-order expansion 
for the matrix element of Z between states v and w of 
an atom with one valence electron: 

(* W \Z\9 V ) =ZW+ZW+--- , (53) 

The first-order matrix element is given by the DHF value 
in the present case 

Z$ = *W (54) 

The second-order expression for the matrix element 
of a one-body operator Z in a Hartree-Fock potential 
is given by 

r7(2) \ * ZamQwmva ^ ^ ZmaQwavm 

z wv = 2^ — + 2^ — - — , (55) 

am £av £mw am €wa £mv 

where e wa = e w + e a . The summation index a ranges over 
states in the closed core, and the summation index m 
ranges over the excited states. The complete third-order 
MBPT expression for the matrix elements of monovalent 
systems was given in Ref. (5lj . The monumental task of 
deriving and evaluating the complete expression for the 
fourth- order matrix elements has been carried out for Na 
in Ref. [H|. 

The polarizabilities are obtained using a sum-over- 
state approach by combining the resulting matrix ele- 
ments and either experimental or theoretical energies. 
The calculations are carried out with a finite basis set, 
resulting in a finite sum in the sum-over-state expres- 
sion that it is equivalent to the inclusion of all bound 
states and the continuum. Third-order MBPT calcula- 
tion of polarizabilities is described in detail, for example, 
in Ref. OH for Yb+. 

The relativistic third-order many-body perturbation 
theory generally gives good results for electric-dipole 
(El) matrix elements of lighter systems in the cases when 
the correlation corrections are not unusually large. For 
example, the third-order value of the Na 3s— 3pi/2 matrix 
element agrees with high-precision experiment to 0.6% 
[37]]. However, the third-order values for the 6s — Gp±/2 
matrix element in Cs and 7s — 7pi/2 matrix element in 



Fr differ from the experimental data by 1.3% and 2%, 
respectively [37| . For some small matrix elements, for 
example 6s — 7p in Cs, third-order perturbation theory 
gives much poorer values. As a result, various methods 
that are equivalent to summing dominant classes of per- 
turbation theory terms to all orders have to be used to 
obtain precision values, in particular when sub-percent 
accuracy is required. 

The relativistic all-order correlation potential method 
that enables efficient treatment of do mina nt core- valence 
correlations was developed in Ref. 125]. It was used 
to study fundamental symmetries in heavy atoms and to 
calculate atomic properties of alkali-me tal atoms and iso- 
electronic ions (see, for example, Refs. [I26lll27l | and ref- 
erences therein). In the correlation potential method for 
monovalent systems, the calculations generally start from 
the relativistic Hartree-Fock method in the V N-1 ap- 
proximation. The correlations are incorporated by means 
of a correlation potential £ defined in such a way that 
its expectation value over a valence electron wave func- 
tion is equal to the RMBPT expression for the correla- 
tion correction to the energy of the electron. Two classes 
of higher-order corrections are generally included in the 
correlation potential: the screening of the Coulomb in- 
teraction between a valence electron and a core electron 
by outer electrons, and hole-particle interactions. Lad- 
der diagrams were included to all orders in Ref. [l28j |. 
The correlation potential is used to build a new set of 
single-electron states for subsequent evaluation of vari- 
ous matrix elements using the random-phase approxima- 
tion. Structural radiation and the normalization correc- 
tions to matrix elements are also incorporated. This ap- 
proach was used to evaluate black-body r adiat ion shifts 
in microwave frequency standards in Refs. |l29l . ll3C| (see 
Section IWC51) . 

Another class of the all-order approaches based on the 
coupled-cluster method is discussed in the next subsec- 
tion. 



F. Coupled-cluster methods 

In the coupled-cluster method, the ex act m any- body 
wave function is represented in the form jl3l| 

I*) =cxp(5)|* (0) ), (56) 

where l^ )) is the lowest-order atomic wave function. 
The operator S for an N-electron atom consists of "clus- 
ter" contributions from one-electron, two-electron, ■ ■ • , 
N-electron excitations of the lowest-order wave function 
|*(°)>: S = Si + S 2 + ■ ■ ■ + S N . In the single-double 
approximation of the coupled-cluster (CCSD) method, 
only single and double excitation terms with Si and S2 
are retained. Coupled-cluster calculations which use a 
relativistic Hamiltonian are identified by a prefix of R, 
e.g. RCCSD. 

The exponential in Eq. (|56|) . when expanded in terms 
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of the n-body excitations S ra , becomes 



form 



= u + s t + s 2 + s 3 



2^ 



S1S2 



I* 



(o)\ 



(57) 

Actual implementations of the coupled-cluster ap- 
proach and subsequent determination of polarizability 
vary significantly with the main source of variation being 
the inclusion of triple excitations or non-linear terms and 
use of different basis sets. These differences account for 
some discrepancies between different coupled-cluster cal- 
culations for the same system. It is common for triple ex- 
citations to be included perturbatively. In this review, all 
coupled-cluster calculations that include triples in some 
way are labelled as CCSDT (or RCCSDT, RLCCSDT) 
calculations with no further distinctions being made. 

We can generally separate coupled-cluster calculations 
of polarizabilities into two groups, but note that details 
of calculations vary between different works. Implemen- 
tations of the CCSDT method in the form typically used 
for the quantum chemistry calculations use gaussian type 
orbital basis sets. Care should be taken to explore the 
dependence of the final results on the choice and size of 
the basis set. The dependence of the dipole polarizabil- 
ity values on the quality of the basis set used has been 
discussed, for example, in Ref. 35]. In those calcula- 
tions, the polarizabilities a re gener ally calculated using 
the finite-field approach [35l l36l Il32 1 . Consequently, such 
CC calculations are not restricted to monovalent systems, 
and RCC calculations of polarizabilit ies of div alent sys- 
tems have been reported in Refs. [H, Il33l Il34| . 

The second type of relativistic coupled-cluster calcu- 
lations is carried out using the linearized variant of the 
coupled-cluster method (referred to as the relativistic all- 
order method in most references), which was first devel- 
oped for a tomic physics calculations and applied to He in 
Ref. 135] . The extension of this met hod to monovalent 
systems was introduced in Ref. [l36j |. We refer to this 
approach as the RLCCSD or RLCCSDT method 
We note that RLCCSDT method includes only valence 
triples using perturbative approach. As noted above, all 
CC calculations that include triples in some way are la- 
belled as CCSDT. The RLCCSDT method uses finite ba- 
sis set of B-splines rather than gaussian orbitals. The B- 
spline basis sets are effectively complete for each partial 
wave, i.e. using a larger basis set will produce negligible 
changes in the results. The partial waves with / = — 6 
are generally used. Third-order perturbation theory is 
used to account for higher partial waves where necessary. 
Very large basis sets are used, typically a total of 500 - 
700 orbitals are included for monovalent systems. There- 
fore, this method avoids the basis set issues generally 
associated with other coupled-cluster calculations. The 
actual algorithm implementation is distinct from stan- 
dard quantum chemistry codes as well. 

In the linearized coupled-cluster approach, all non- 
linear terms are omitted and the wave function takes the 



|tf ) = {1 + Si + S 2 + S 3 + ■ ■ ■ + S N } |*( 0) ) . (58) 

The inclusion of the nonlinear terms withi n the frame- 
work of this method is described in Ref. [l37j | . Restricting 
the sum in Eq. (155)) to single, double, and valence triple 
excitations yields the expansion for the wave function of 
a monovalent atom in state v: 



1 + E 

Pma^fn^a / Pm.nab^ rn Ci n Clb^a'^ 

ma mnab 

^ ^ PmvCl> m Clv ~t~ ^ ^ Prnnva^rn^ii^a-^v 
m/t) mna 



- E 

6 ^ 



P m n r v a b Q> m & n Oj \ CI b CI a CI v 



mnrab 



|*<°>>, (59) 



where the indices m, n, and r range over all possible vir- 
tual states while indices a and b range over all occupied 
core states. The quantities p ma , p m v are single-excitation 
coefficients for core and valence electrons and p mna b and 
Pmnva are double-excitation coefficients for core and va- 
lence electrons, respectively, p m nrvab are the triple va- 
lence excitation coefficients. For the monovalent systems, 
U is generally taken to be the frozen- core V^- 1 potential, 
U = V DF . 

We refer to results obtained with this approach as RL- 
CCSDT, indicating inclusion of single, double, and par- 
tial triple excitations. The triple excitations are generally 
included perturbatively. Strong cancellations between 
groups of smaller terms, for example non-linear terms 
and certa in triple excitation terms have been found in 
Ref. [138| . As a result, additional inclusion of certain 
classes of terms may not necessarily lead to more accu- 
rate values. 

The matrix elements for any one-body operator Z 
given in second-quantized form by Eq. (|52j) are obtained 
within the framework of the linearized coupled-cluster 
method as 



y/(* v \V v )(* w \V* 



(60) 



where \^ v ) and \^ w ) are given by the expansion (|59|) . In 
the SD approximation, the resulting expression for the 
numerator of Eq. (l60l) consists of the sum of the DHF 
matrix element z wv and 20 other terms that are linea r 
or quadratic functions of the excitation coefficients ( 1 36l j . 
The main advantage of this method is its general applica- 
bility to calculation of many atomic properties of ground 
and excited states: energies, electric and magnetic mul- 
tipole matrix elements and other transition properties 
such as oscillator strengths and lifetimes, A and B hy- 
pcrfmc constants, dipole and quadrupole polarizabilities, 
parity-nonconserving matrix elements, electron electric- 
dipole-moment (EDM) enhancement factors, C3 and Cq 
coefficients, etc. 
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The all-order method yields results for the properties of 
alkali atoms [31( in excellent agreement with experiment. 
The application of this method to the calculation of al- 
kali polarizabilities (using a sum-over-state approach) is 
described in detail in Refs. [H-fHHl]. 

In its present form described above, the RLCCSDT 
method is only applicable to the calculation of polariz- 
abilities of monovalent systems. The work on combining 
the RLCCSDT approach with the CI method to cre ate a 
method that is more general is currently in progress [l39j | 
and is described in Section IIVHI 



the two-body part of core-valence interaction that rep- 
resents screening of the Coulomb interaction by valence 
electrons; H 2 — > H 2 + £2 - Th e quantity E 2 is calculated 
in second-order MBPT [lOOj]. The CI method is then 
used with the modified i? c ff to obtain improved energies 
and wave functions. 

The polarizabilities are determined using the direct ap- 
proach (in the valence sector) by solving the inhomoge- 
neous equation in the valence space, approximated from 
Eq. ([6]) . For state v with total angular momentum J and 
projection M, the corresponding equation is written as 



G. Combined CI and many-body perturbation 
theory 

Precise calculations for atoms with several valence elec- 
trons require an accurate treatment of valence-valence 
correlations. While finite-order MBPT is a powerful tech- 
nique for atomic systems with weakly interacting con- 
figurations, its accuracy can be limited when the wave 
function has a number of strongly interacting configura- 
tions. One example occurs for the alkaline-earth atoms 
where there is strong mixing between the ns 2 and np 2 
configurations of 1 S symmetry. For such systems, an ap- 
proach combining both aspects has been developed by 
Dzuba et al. [lOOj ] and later applied to the calculati on o f 
atom ic properties of many other systems [H, H3, Il40l — 
Il43l |. This composite approach to the calculation of 
atomic structure is often abbreviated as CI+MBPT (we 
use RCI+MBPT designations in this review to indicate 
that the method is relativistic). 

For systems with more than one valence electron, the 
precision of the CI method is drastically limited by the 
sheer number of the configurations that should be in- 
cluded. As a result, the core-core and core- valence corre- 
lations might only receive a limited treatment, which can 
lead to a significant loss of accuracy. The RCI+MBPT 
approach provides a complete treatment of core corre- 
lations to a limited order of perturbation theory. The 
RCI+MBPT approach uses perturbation theory to con- 
struct an effective core Hamiltonian, and then a CI cal- 
culation is performed to generate the valence wave func- 
tions. 

The no-pair Hamiltonian given by Eqs. (f4"5j) and ([4"6")) 
separates into a sum of the one-body and two-body in- 
teractions, 



(E v - flee)|*(«,Af')> = Arfr|*o(», J,M)). 



(62) 



H = H 1 +H 2 



(61) 



where H 2 contains the Coulomb (or Coulomb + Breit) 
matrix elements Vijki- In the RCI+MBPT approach, the 
one-body term Hi is modified to include a correlation 
potential Ei that accounts for part of the core-valence 
correlations, Hi — > Hi + Ex. Either the second-order 

(2) 

expression for or all-order ch ains of such terms can 
be used (see, for example, Ref. jlOOl ]). The two-body 
Coulomb interaction term in H 2 is modified by including 



The wave function &(v, M') is composed of parts that 
have angular momenta of J' = J, J±l. This then permits 
the scalar and tensor polarizability of the state \v,J, M) 
to be determined [53j . 

The construction of H e g was described in the preceding 
paragraphs. The effective dipole operator D e g includes 
random phase approximation (RPA) correcti ons a nd sev- 



144 . Non- 



eral smaller MBPT corrections described in 
RPA corrections may be neglected in some cases [53|. 
There are several variants of the RCI+MBPT method 
that differ by the corrections included in the effective op- 
erators -ffeff and -D e ff, the functions used for the basis 
sets, and versions of the CI code. In some implemen- 
tations of the RCI+MBPT, the strength of the effective 
Hamiltonian is rescaled to improve agreement with bind- 
ing energies. However, this procedure may not necessar- 
ily improve the values of polarizabilities. 

The contributions from the dominant transitions may 
be separated and replaced by more accurate experimental 
matrix elements when appropriate. Such a procedure is 
discussed in detail in Ref. [14lj |. This hybrid RCI+MBPT 
approach [rJ Il45| has been used to obtain present 
recommended values for the polarizabilities of the ns 2 
and nsnp 3 Pq states of Mg, Ca, Sr, Hg, and Yb needed 
to evaluate the blackbody radiation shifts of the relevant 
optical frequency standards. 



H. Combined CI and all-order method 

The RCI+MBPT approach described in the previous 
section includes only a limited number of the core- valence 
excitation terms (mostly in second order) and deterio- 
rates in accuracy for heavier, more complicated systems. 
The linearized coupled-cluster approach described in Sec- 
tion IIVFI is designed to treat core-core and core- valence 
correlations with high accuracy. As noted above, it is re- 
stricted in its present form to the calculation of properties 
of monovalent systems. Direct extension of this method 
to even divalent systems faces two major problems. 

First, use of the Rayleigh-Schrodinger RMBPT for 
heavy systems with more than one valence electron leads 
to a non-symmetric effective Hamiltonian and to the 
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problem of "intruder states" [1461 ] . Second, the com- 
plexity of the all-order formalism for matrix elements in- 
creases rapidly with the number of valence electrons. The 
direct extensions of the all-order approach to more com- 
plicated systems is impractical. For example, the expres- 
sion for all-order matrix elements in divalent systems con- 
tains several hundred terms instead of the twenty terms 
in the corresponding monovalent expression. However, 
combining the linearized coupled-cluster approach (also 
referred to as the all-order method) with CI method elim- 
inates many of these difficulties. This m etho d (referred 
to as CI+all-order) was developed in Rcf. 139] and tested 
on the calculation of energy levels of Mg, Ca, Sr, Zn, Cd, 
Ba, and Hg. The prefix R is used to indicate the use of 
the relativistic Hamiltonian. 

In the RCI+all-order approach, the effective Hamilto- 
nian is constructed using fully converged all-order exci- 
tations Coefficients Pma> Pmriabi Pmv i Pmnvai and Pmnvw 

(see section IIVFI for designations) . The p mnV w coeffi- 
cients do not arise in the monovalent all-order method, 
but are straightforwardly obtained from the above core 
and core-valence coefficients. As a result, the core-core 
and core- valence sectors of the correlation corrections for 
systems with few valence electrons are treated with the 
same accuracy as in the all-order approach for the mono- 
valent systems. The CI method is used to treat valence- 
valence correlations and to evaluate matrix elements and 
polarizabilities. 

The RCI+all-order method employs a variant of the 
Brillouin-Wigner many-body perturbation theory, rather 
than Rayleigh-Schrodinger perturbation theory. In the 
Brillouin-Wigner variant of MBPT, the effective Hamil- 
tonian is symmetric and accidentally small denominators 
do not arise [l39j ] . Comparisons of the RCI+MBPT and 
RCI+all-order binding energies for the ground and ex- 
cited states of a number of two-electron systems reveal 
that the RCI+all-order energies are usually more accu- 
rate by at least a factor of three [l39j ] . 

The preliminary calculations of polarizabilities values 
in Ca and Sr indicate better agreement of the RCI+all- 
order ab initio results with recommended values from 
Ref. [13| in comparison with the RCI+MBPT approach. 



trie constant, RI - refractive index, SA - spectral analysis, 
RRPA - relativistic random phase approximation, MBPT 
- many-body perturbation theory, (R)CCSDT - (rela- 
tivistic) coupled-cluster method. If any triple excitations 
are included, CCSDT abbreviation is used for coupled- 
cluster calculations, single-double coupled cluster calcu- 
lations are labelled (R)CCSD. The RCCR12 calculation 



1511 ] is a CCSDT calculation which allows for explic- 



itly correlated electron pairs. The pseudo-natural orbital 
coupled electron pair approximation (PNO-CEPA) [104| 
can be regarded as precursor of modern CCSD type mod- 
els. We first discuss the general trends of values for the 
noble gases as a whole, and then consider He in more 
detail separately. 

The most precise calculations of the noble gas polariz- 
abilities (apart from helium) have mostly been obtained 
with coupled-cluster type calculations. As we noted in 
the previous sections, particular care has to be taken to 
ensure that the basis set used in CC calculations is of 
sufficiently high quality to obtain accurate values. One 
curious aspect about the noble gases is their insensitivity 
to relativistic effects. The relativistic correction to ao is 
less th an 1% for Ne, Ar, and Kr and is only about 2% 
for Xe [III. 



One notable feature of Table IIVI is the good agree- 
ment of the RRPA [50| with the much more elaborate 
coupled-cluster and Hylleraas basis function calculations 
and experimental data. The difference between RRPA 
values and other calculations/experimental value for neu- 
tral systems ranges from 10% for Ne to 1.6% for Kr (4% 
for He). The RRPA values [5(| improve significantly for 
the singly ionized systems and differ from other values by 
5% for Na + and only 0.4% for Rb + . The discrepancies 
are reduced further for doubly ionized systems owing to 
the decrease in the relative contribution of the correla- 
tion corrections beyond RRPA. Core polarizabilities for 
the alkali and alkaline-earth atoms are important for the 
construction of CICP type models of these atoms. In ad- 
dition, the RRPA calculations of the core polarizabilities 
are embedded into many calculations of the polarizabil- 
ities of alkali and alkaline-earth ions (see, for example, 
Refs. [llllllET 



BENCHMARK COMPARISONS OF THEORY 
AND EXPERIMENT 

A. Noble gases and isoelectronic ions 



Theoretical M. M. Hoi fTT6l fTl7l .^147L /l48l [l50l [l5Tl 
158 . l280j ] and experimental "p. M I1T2T I152M155L I159I- 
1671 ] values for the ground state polarizabilities of the 
noble gases and isoelectronic ions are listed in Table IIVI 
References are given in square brackets. The reference 
is given at the end of the row when all data in this row 
come from the same work. Otherwise, the references are 
listed together with the particular value. The following 
method abbreviations are used in the table: DC - dielec- 



1. Helium 

The helium atom is of particular interest since it al- 
lows for the most precise calculations and benchmark 
tests of theory and experiment. Within the frame- 
work of the non-relativistic Schrodinger equation with 
infinite-nuclear-mass Hamiltonian, the He polarizability 
value obtained using a modified version of the gener alized 
Hylleraas basis set [281] is 1.383192174455(1) a.u. [Ill], 
achieving accuracy of 13 significant di gits. This value is 
in agreement with 1996 calculation of |117| . 

The finite mass effects increase the polarizability by 
about 0.00062 a.u., with the mass polarization effect ac- 
counting for 0.000049 a.u. resulting in the 4 He nonrel- 
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TABLE IV: Ground state polarizabilities Qo (in atomic units) of noble gases and isoelectronic ions. Uncertainties in the last 
digits are given in parentheses. References are given in square brackets. Method abbreviations: DC - dielectric constant, RI - 
refractive index, SA - spectral analysis, RRPA - relativistic random-phase approximation, (R)CCSDT - (relativistic) coupled- 
cluster calculations. The RCCR12 calculation is a CCSDT calculation which allows for explicitly correlated electron pairs. "See 
text for further discussion of He polarizability calculations, b Finite mass Hylleraas alculation incorporating relativistic effects 
from an RCI calculation as an additive correction, c PNO-CEPA (pseudo-natural orbital coupled electron pair approximation). 



He 


Ne 


Ar 


Kr 


Xe 


Rn 


Method [Ref.] 


1.322 


2.38 


10.77 


16.47 


26.97 




Th. RRPA [50] 


1.383763 


2.6648 


11.084 








Th. CCSDT [147] 




2.697 


11.22 


16.80 


27.06 


33.18 


Th. RCCSDT [148] 




2.665 [1491 


11.085(6) [361 








Th. CCSDT 




2.6557 


11.062 


17.214 


28.223 




Th. MBPT [150] 




2.668(6) [151] 










Th. RCCR12 


1.38376079(23) a ' 6 [119] 












Th. 


1.383223(67) [152. 1531 


2.670(3) [1541 


11.081(5) [154] 


16.766(8) [1541 






Expt. DC 


1.3838 


2.6680 


11.091 


16.740 


27.340 




Expt. RI [62] 


1.384 


2.663 


11.080 


16.734 


27.292 




Expt. RI [61] 


1.383759(13) [63] 




11.083(2) [1551 








Expt. RI 


Li+ 


Na+ 


K+ 


Rb+ 


Cs+ 


Fr+ 




0.192486 6 [156. 1571 


0.9947 c [1041 


5.354 c [104] 








Th. 


0.1894 


0.9457 


5.457 


9.076 


15.81 




Th. RRPA [50] 




1.00(4) 


5.52(4) 


9.11(4) 


15.8(1) 


20.4(2) 


Th. RCCSDT [158] 


0.1883(20) [1591 


0.978(10) [160] 


5.47(5) [160] 


9.0 [161] 


15.544(30) [162] 




Expt. SA 




1.0015(15) [163] 






15.759 [164] 




Expt. SA 




0.9980(33) [1651 






15.644(5) [112. 1661 




Expt. SA 


Be 2+ 


Mg 2+ 


Ca 2 + 


Sr 2+ 


Ba 2+ 


Ra 2+ 




0.05182 


0.4698 


3.254 


5.813 


10.61 




Th. RRPA [50] 


0.052264 6 [156. 1571 


0.4814 c [1041 


3.161 c [1041 








Th. 






3.262 


5.792 


10.491 


13.361 


Th. RCCSDT [35] 




0.489(5) [1601 


3.26(3) [1601 








Expt. SA 




0.486(7) [1671 










Expt. SA 



ativistic value of 1.38380999 a.u. [116|, lllg l282j . The 
a 2 r e lativ i stic corrections contribute -0.00008035(2) a.u. 
[HI EH |Hi|. The a 3 QED corrections with excep- 
tion of the terms containing electric-field deriv ative of 
the Bethe logarithm were calculated in Ref. |ll6l ] to give 
0.0000305 a.u.. These latter terms were calculated in 
[l 19j j . together with the estimates of the a 4 , a 2 m e /Mn c , 
and a 3 m e /MHe, yielding the final value of 4 He polariz- 
ability of 1.38376079(23) a.u. listed in Table [TV] 

A non-relativistic coupled-cluster c alcu lation in the in- 
finite mass limit carried out in Ref. |l47l | provides a de- 
tailed study of the dependence of the CCSDT results on 
the choice of the basis set and tests of basis set conver- 
gence. The values obtained with different uncontracted, 
even-tempered basis sets varied in the fifth significant 
digit. Their final value of a(°°He) = 1.383763 a.u. dif- 
fers from the exa c t no n-relativistic Hylleraas value of 
1.383192 a.u. [HI, [H3 at the same level. 

A microwave cavity was recently used to measure 
the refractive index of helium giving a polarizability of 
1.383759(13) a.u. (63|. The best experiment has an un- 



certainty of about 10 p pm a nd is in accord with the most 
accurate theory value |l!9j |. Availability of such precise 
theoretical and experimental values of He polarizability 
allows for accurate determinations of the thermodynamic 
temperature and may lead to a more accurate value of the 
Boltzmann constant 63] . This application is discussed in 
more detail in Section IVHEl 



B. Monovalent systems 



T he theore tical 

Ell EMJrn] 



7l} (l78l liia 




M, EM EH EM 

Il88l4l93l 12831 a nd experimen- 
tal [67J, 17114751 I90L Ti811 Il83l4l86l ll94H199l [284] values of 
static scalar (ao) and tensor (012) polarizabilities of al- 
kali atoms and scalar static polarizabilities of singly ion- 
ized monovalent ions are compared in Table |V] and I VII 
The same designations are used as in the noble gas ta- 
ble. The following additional method abbreviations are 
used: EH - E-H balance or beam-deflection, sum-rule - 
hybrid /-sum rules with experimental data for primary 
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TABLE V: Ground and npj excited state polarizabilities (in a.u.) of alkali atoms. Scalar (qo) and tensor (02) polarizabilities 
are given for the np 3 /2 states. Static polarizabilities for the npi/2 and np 3 /2 states are the same for the non-relativistic 
Hylleraas and CICP calculations. Uncertainties in the last digits are given in parentheses. References are given in square 
brackets. Method abbreviations: EH - E-H balance or beam-deflection, sum-rule - hybrid /-sum rules with experimental data 
for primary contribution, SA - spectral analysis, CI - configuration interaction, CICP - CI calculations with a semi-empirical 
core potential, MBPT - many-body perturbation theory, RLCCSDT - linearized CCSD method with partial triple contributions. 
All values in the sum-rule row explicitly include a core polarizability. a Non-relativistic Hylleraas calculation for °°Li, b Hylleraas 
calculations for 7 Li that includes estimate of relativistic effects, C CI, d Hybrid-RLCCSD data for the alkali ground states from 
[56| are listed as recommended "sum-rule" data, e interferometry, ^interferometry ratio, 9 cold atom velocity change experiments. 



Li 


Na 


K 


Rb 


Cs 


Fr 


Method 




qo 2s 


3s 


4s 


5s 


6s 


7s 






164.112(l) a [118] 164. 50 c [168] 






398.2(9)" [169] 




Th. 




164.11(3) 1 170 J 












Th. Hyl. 




164.21 11711 


162.8 I105J 
165.50 11721 


290.0 [105] 
301.28 [173] 


315.7 [105] 






Th. CICP 
Th. CCSD 




1 f?o 7/1 r-i 701 

IDo. / 4 ^^^] 




291.12 [1321 


316.17 [1321 


396.02 [1321 


315.23 [132] Th. RCCSDT 






IDo.U 


289.1 


316.4 


401.5 


315.1 


Th. RLCCSD [56] 




no M 7K1 
164.08 llYOl 




289.3 [1761 




398.4(7) [177] 


313.7 [1781 


Th. RLCCSDT 




1 C A ( O \ 

164(oj 


159(3) 


293(6) 


319(6) 


402(8) 




Expt. EH [67J 




164.2(1. l) e [71] 


162.7(8) e [72J 


290.8(1.4) / [73J 


318.8(1.4) / [73] 401.0(6) 9 [74J 




Expt. 






162.6(3) 


290.2(8) 


318.6(6) 


399.9(1.9) 


317.8(2.4) 


Sum-rule 1 * [56] 




Qo 2pi /2 


3P1/2 


#1/2 


5^1/2 


6f>l/2 








126.9458(3) [118] 












Th. a Hyl. 




126.95 [171] 


360.7 [1791 


615.3 [1791 


854.4 [179] 






Th. CICP 




126.980 [175] 




604.1 [1761 


805(31) [30J 


1338(54) [1771 




Th. RLCCSDT 






359.7 


605 


807 






Th. RCI+MBPT 


180] 


Qo 2p 3/2 


3P3/2 


4P3/2 


5^3/2 


6P3/2 








126.995 [175] 




614.1 [1761 




1648(58) [1771 




Th. RLCCSDT 






361.4 


616 


870 






Th. RCI+MBPT 


180] 


Q2 2p 3/2 


3P3/2 


4f>3/2 


5f3/2 


6f>3/2 








1.6214(3) [118] 












Th. Hyl. 




1.6627 [171] 


-87.89 [1791 


-107.9 [1791 


-160.5 [1791 






Th. CICP 






-88.0 


-111 


-171 






Th. RCI+MBPT 


180] 


1.59 [175] 




-107.9 [176] 




-261(13) [177] 




Th. RLCCSDT 




1.64(4) [181] 


-88.3(4) [182] 
-113(16) [1841 


-107(2) [183] 
-110.9(2.8) [1851 


-163(3) [183] 


-261(8) [90] 
-262.4(1.5) [1861 




Expt. 
Expt. 





contribution, RESIS - resonant excitation Stark ioniza- 
tion spectroscopy, RLCCSDT - linearized CCSD method 
with partial triple contributions included. First, some 
general remarks are made for monovalent systems, and 
then Li, Na, Mg + and Cs are considered in more detail. 

The comparatively simple electronic structure of these 
atoms render them amenable to accurate calculation by 
the coupled-cluster and CICP methods. The sum-rule 
polarizabilities |5a ] come from a hybrid calculation that 
use the RLCCSD calculation as a template. However, 
the matrix element for the resonance transition has been 
replaced by high accuracy experimental matrix elements 
compiled in [3l|. The ab initio RLCCSD values are in 
excellent agreement (better than 1%) with these hybrid 
recommended values. The semi-empirical CICP calcula- 
tions reveal a similar level of accuracy, although there 



has been some degradation in accuracy for the heavier 
Rb system. The CI calculations with a semi-empirical 
core potential (CICP) are in excellent agreement with 
RLCCSDT calculations and experiment for lighter sys- 
tems. The non-relativistic CICP cannot be expected to 
be particularly accurate for states with significant spin- 
orbit splitting, e.g. the npj states of Rb. The best that 
can be expected is that the CICP calculation will do a 
reasonable job of reproducing the statistically weighted 
npj average polarizability. 

The results of the coupled-cluster calculations can 
be sensitive to particular contributions that are in- 
cluded, owing to cancellations of various terms (for ex- 
ample, some triple excitations beyond perturbative treat- 
ment may partially cancel with non-linear single-double 
terms), leading to some differences between different 
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TABLE VI: Ground state polarizabilities (in a.u.) of alkali-like ions. Uncertainties in the last digits are given in parentheses. 
References are given in square brackets. Method abbreviations: SA - spectral analysis, RESIS - resonant excitation Stark 
ionization spectroscopy, °non-relativistic Hylleraas calculation for °°Be + , ^Hylleraas calc ulat ions for 9 Be + that includes estimate 
of relativistic effects, c /-sum rule for valence polarizability with core-polarzation from [1871 ) added. 



Be+ 


Mg+ 


Ca+ 


Sr+ 


Ba+ 




Ra+ 


Method 


2s 


3s 


4s 


5s 


6s 




7s 




24.4966(l) a [188] 














Th. Hyl. 


24.489(4) b [170] 














Th. Hyl. 


24.495 [189] 


35.66 [168] 












Th. CI 


24.493 [1881 


34.99 [84J 


75.49 [190] 


89.9 [1911 








Th. CICP 




35.05 [84J 


76.1(1.1) [1921 


91.3(9) [1931 


124.15 [96 




106.5 [1781 


Th. RLCCSD 






75.88 


91.10 


123.07 




105.37 


Th. RCCSDT [35] 




33.80(50) [194] 


75.3 [195] 




125.5(1.0) 


[196] 




Expt. SA 




35.04(3) [84J 






124.30(16) 


[197] 




Expt. RESIS 




35.00(5) [198J 






123.88(5) 


199] 




Expt. RESIS 




35.10 [200] 


74.11 [200] 










c /-sum rule 



TABLE VII: Polarizability differences Qo(npj) — Qo(ns) (in a.u.) of the alkali atoms derived from Stark shift measurements. 
Values are negative when the npj state polarizability is smaller than the ground state polarizability. Stark shifts for the npi/ 2 
and states are the same for the non-relativistic Hylleraas and CICP methods. Uncertainties in the last digits are given 

in parentheses. References are given in square brackets. The experimental values and Hylleraas calculations [170| ] are those 
reported for 7 Li, the CICP and RLCCSDT values are for °°Li. 



7 Li 


Na 


K 


Rb 


Cs 


Method 


2s-2pi /2 


3s-3pi /2 


4s-4p 1/2 


5s-5p 1/2 


6s-6pi /2 




-37.14(3) [1701 










Th. Hylleraas 


-37.26 [105. 1711 


197.9 [105, 1791 


325.3 [105. 1791 






Th. CICP 


-37.104 [175] 


196.7 [56. 1801 


314.8 [1761 


488(4) 6 [30. 1801 


940(55) 6 [177J 


Th. RLCCSDT 


-37.146(17) [91] 




316.68(4) [201] 


491.52(6) [2011 


926.08(12) [92] 


Expt. 


-37.11(33) [1811 


196.86(45) [2021 


315(3) [1831 






Expt. 


2s-2p 3/2 


3s-3p 3/2 


4s-4p 3/2 


5s-5p 3/2 


6s-6p 3/2 




-37.089 MBPT [175] 


198.4 [56, 180] 


324.8 [176] 


554 [56, L80J 


1250(59) [177] 


Th. RLCCSDT 


-37.30(42) 11811 


198.0(6) 11821 


322.3(3.2) [1831 


538.5(3.2) [1831 


1240.2(2.4) [186] 


Expt. 










1264(13) [90] 


Expt. 



coupled-cluster calculations 285] . The properties involv- 
ing nd states (i.e. np polarizabilities) are also sensitive 
to the number of partial of waves included in the basis 
sets. Omission or inadequate inclusion of partial waves 
with I > 3 may lead to poor results for matrix elements 
involving nd states, and, subsequently, relevant excited- 
state polarizabilities. 

Some of the most stringent tests of polarizability cal- 
culations of monovalent systems come from Stark shift 
measurements of alkali resonance transitions. Therefore, 
it is useful to compare the experimental values for the po- 
larizability difference obtained from Stark shift measure- 
ment directly with theoretical predictions in these cases. 
Scalar polarizability differences ao(npj)—ao(ns) (in a.u.) 
of the alkali transitions derived from Stark shift mea- 
surements are co mpared with theoretical values in Ta - 
ble ryiii f3ol [56l M M. H05L [1701 [ml EMOxE EM E13- 

ll83jl8aT2CnTl202| . For the elements heavier that Li, the 



finite mass effects are smaller than the uncertainty of the 
calculation. 



The tensor polarizability of an open shell atom can be 
extracted from the difference in polarizabilities between 
the different magn etic sub- l evels. The scalar and ten- 
sor polarizabilities [ill Ell EaMll HqUqI, SqI [208L 
l209j ] of some low lying excited states of Li, Na, K, Rb, 
Ca+ and Sr + are listed in Table IVIIII There is a paucity 
of experimental data for excited states, even for well- 
studied alkali atoms. The polarizabilities of the nd^/2 
states of Ca + and Sr + are given owing to their impor- 
tance for evaluation of the black-body radiation shifts. 
Some older and less accurate Stark shifts and tensor po- 
larizabilities are omitted from these Tables. 
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TABLE VIII: Excited state scalar cto and tensor 0:2 polar- 
izabilities (in a.u.) of monovalent systems. All experimental 
values are derived from Stark shift experiments and the polar- 
izability of the lower state is added to the Stark shift to get the 
upper state polarizability. Uncertainties in the last digits are 
given in parentheses. References are given in square brackets. 
a Hylleraas basis functions, °CICP, C RLCCSDT, d CA. Polar- 
izabilities marked with an asterisk (*) were not published, but 
obtained from the matrix elements of [1791 ]. 



Atom 


State 




Expt. 




Theory 


Li 


3d 3 /2 


«() 


-15082(60) 


[2031 


-14928" 


[118] 












-15044 6 


[171] 






a.i 


11626(68) 


12031 


11409 a 


[1 1 81 












11490 6 


[1711 




3rfr In 


Q() 


-15159(32) 


[2031 


-14928° 


[118] 












-15044 6 


[1711 






Ct2 


16308(52) 


[2031 


16298° 


[118] 












16414 6 


[171] 


Na 




OCQ 


21000(1200) 


[2041 


21780* 6 


[1791 




O / z 


Q() 


624000(7000) 


[2041 


633800* 6 


[179] 






Ct2 


-154700(2800) 


[2041 


-148700* 6 


[179] 




4rf 5 /2 


cto 


627000(5000) 


[2041 










OL1 


-213800(2000) 


[204] 


-212400* 6 


[179] 


K 


5^3/2 


Cto 






7118* 6 


[1791 






Q •> 


-1057(161) 


[2051 


-1019* 6 


[1791 


Ca+ 


3rf 5 /2 


cto 






32. 73 6 


[190] 












32.0(1. l) c 


[192] 




3rf 5 /2 


«2 






-25.20 6 


[190] 












-24.5(4) c 


[192] 


Rb 


6P3/2 


a>2 


-2090(80) 


[2061 


-2040 d 


[2071 




6d 3 / 2 


a 2 


-42.2(28) 


[208] 


-559 d 


[2071 




6d 5 / 2 


a>2 


3780(200) 


[2081 


3450 d 


[2071 




7P3/2 




-12900(800) 


[209] 


-12500 d 


[207] 


Sr+ 


4rf 5 /2 


Cto 






61.77° 


[191J 












62.0(5) c 


[193] 




4rf 5 /2 


a 2 






-47.20° 


[191J 












-47.7(3) c 


[193] 



Lithium 

The lithium polarizability could assume a pivotal role 
in polarizability metrology if a multi-species interferom- 
eter can be constructed that is capable of measuring the 
ratio of the polarizability of other atoms to that of Li 
to a relative accuracy of 10 -4 [7(|. In this case, a mea- 
surement of such ratios in conjunction with a definitive 
calculation of the Li cvo could lead to new accuracy bench- 
marks for the polarizabilities of a number of elements. 

Correlated basis calculations are possible for lithium 
since it only has three electrons. Consequently it has 



TABLE IX: Selected theoretical and experimental ground 
state polarizabilities cto (in a.u.) of sodium atom. Uncer- 
tainties in the last digits are given in parentheses. References 
are given in square brackets. HF - Hartree-Fock, PNO-CEPA 

- pseudonatural orbital configuration expansion, CICP - CI 
calculations with a semi-empirical core potential, RLCCSDT 

- linearized CCSD method with additional partial triple con- 
tributions included, EH - E-H balance or beam-deflection. 



Method 


Year 


Value 




Theory 




HF [9101 


1964 


183 


HF [1041 


1984 


189.2 


PNO-CEPA [2111 


1976 


165.02 


PTPP [9191 


1 Q7Q 


1 fi9 fi 


CICP [1041 


1984 


162.4 


CICP [1051 


2003 


162.8 


CI [1681 


2007 


164.50 


RLCCSD [56] 


1999 


163.0 


RCnSDT [1 731 


1999 


164.89 


CCSDT [213] 


2001 


165.06 


RCCSDT [214] 


2003 


166.3 


RCCSDT [172] 


2004 


165.5 


CCSDT [174] 


2005 


162.88(60) 




Experiment 




/-sum [215] 


1959 


166 


EH [66] 


1974 


165(11) 


EH [67] 


1974 


159(3) 


Interferometry [72J 


1995 


162.7(8) 


Hybrid /-sum [56] 


1999 


162.6(3) 


Interferometry [731 


2010 


162.7(1.3) 



been poss ible to ca lculate the polarizability to very high 
precision [l!7l Ill8j | . The uncertainty in the experimental 
value of the polarizability 164.2(11) a.u. [7l[ spans all of 
the theoretical results reported in Table [V] 

The most recent Hyl leraa s calculation gave ao — 
164.112(1) a.u. for °°Li [HI- Including finite mass ef- 
fects gave ao = 164.161(1) a.u. for 7 Li. An approxi- 
mate treatment of relativistic effects gave a recommended 
value of 164.11(3) a.u. [l70j |. Hylleraas polarizabilities 
could also serve as benchmarks for coupled-cluster type 
calculations which can be applied to atoms heavier than 
lithium. 

The most stringent test of Li polarizability calculations 
is presently the Stark shift measurement of the 2s-2pi/2 
transitions by Hunter et al. (9l| . which gave a polar- 
izability difference of —37.14(2) a.u. The current the- 
oretical benchmark is the recent Hylleraas c alculation s 
that include finite mass and relativistic effects [1181 . 1 1 TOj j . 
The 7 Li H ylleraas polarizability difference of —37.14(4) 
a.u. [l70j | is in excellent agreement with the experimen- 
tal polarizability difference [9l|. The RLCCSDT value of 
—37.104 is within 2 standard deviations of the Hunter ex- 
periment while the CICP value is 4 standard deviations 
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TABLE X: Excited state scalar ao and tensor ai polarizabilities (in multiples of 1000 a.u.) of the Cs atom. Uncertainties 
in the last digits are given in parentheses. References are given in square br acket s. Expe rimental values: der ived fro m th e 
Ref. [p 7 s - 6s Stark shift measurement and the 6s result from c Ref. [21(J, d Refs. HO,!!!!], e Ref. [H3], / Ref. [220( ]. 
9 Ref. 221], h Ref. [222t ] . 'Ref. }223| , CA - Coulomb approximation, RLCCSDT - relativistic linearized coupled-cluster method 
with single, double, and partial triple excitations. 



ao 


7s 


8s 


9s 


10s 


lis 


12s 


Ref. 






6.238(41) 


38.27(28) 


153.7(1.0) 


478(3) 


1246(8) 


2866(30) 


Th. RLCCSDT 


[177] 




6.14 


37.9 


153 


475 


1240 


2840 


Th. CA [2181 






6.238(6/ 


38.06(25) c 




478.5(1.1/ 


1245(1/ 


2867(2) d 


Expt. 




Q'O 


7pl/2 


8P1/2 


9^1/2 


10pi/2 












29.9(7) 


223(2) 


1021(7) 


3499(19) 






Th. RLCCSDT 


[177] 




29.4 


221 


1020 


3490 






Th. CA [2181 






29.6(6) 












Expt. [224J 






7P3/2 


8^3/2 


9P3/2 


10p 3 /2 










ao 


37.5(8) 


284(3) 


1312(7) 


4522(19) 






Th. RLCCSDT 


[177] 




36.9 


282 


1310 


4510 






Th. CA [2181 






37.9(8) 












Expt. [219] 




a 2 


-4.41(17) 


-30.6(6) 


-135(2) 


-451(5) 






Th. RLCCSDT 


[177] 




-4.28 


-30.2 


-134 


-449 






Th. CA [2181 






-4.43(12/ 


-30.5(1.2) / 










Expt. 






-4.33(17) 












Expt. [2061 






-4.00(8) 












Expt. [224] 






5d 3 / 2 


6d 3 / 2 


7rf 3 / 2 


8^3/ 2 


9rf 3/2 


10d 3 /2 






ao 


-0.352(69) 


-5.68(45) 


-66.7(1.7) 


-369(5) 


-1402(13) 


-4234(32) 


Th. RLCCSDT 


[177, 225] 




-0.418 


-5.32 


-65.2 


-366 


-1400 


-4220 


Th. CA [2181 










-60(8/ 




-1450(120) / 


-4185(4/ 


Expt. 




a 2 


0.370(28) 


8.77(36) 


71.1(1.2) 


339(4) 


1189(10)' 


3416(26) 


Th. RLCCSDT 


[177, 225] 




0.380 


8.62 


70.4 


336 


1190 


3410 


Th. CA [2181 










74.5(2.0)' 


332(16/ 


1183(35)' 


3401(4) h 


Expt. 






5d 5/2 


6d 5 / 2 


74/2 


8d 5 /2 


9d 5 /2 


10d 5 /2 






Q'O 


-0.453(70) 


-8.37(55) 


-88.8(2.0) 


-475(5) 


-1777(14) 


-5316(38) 


Th. RLCCSDT 


[177, 225] 




-0.518 


-7.95 


-87.1 


-472 


-1770 


-5300 


Th. CA [2181 










-76(8/ 




-2050(100/ 


-5303(8)'' 


Expt. 




a 2 


0.691(40) 


17.33(50) 


142(2) 


678(5) 


2386(13) 


6869(34) 


Th. RLCCSDT 


[177, 225] 




0.704 


17.00 


140 


675 


2380 


6850 


Th. CA [2181 










129(4) 9 


73 1(40) f 


2650(140) / 


6815(20/ 


Expt. 
















7110(360) 


Expt. [220J 





too large. Table IVTil shows that the Stark shift data offer 
the most precise information with which to discriminate 
between various theoretical calculations. 



2. Sodium 

A chro nolog ical list detail ing selected values [H, |H, 

of the sodium 

ground state polarizability is presented in Table IIX1 The 



theory values are also sorted by the type of calculation. 
The 3s — > 3p resonant transition accounts for 98.8% of 
the polarizability. 

The most notable feature of this table is the excel- 
lent agreement of the semi-empirical CICP type calcula- 
tions with the recent high-precision experimental values 
of 162.6(3) a.u. ^1 and 162.7(8) a.u. All three 

calculations 

[IMlloElIIl, performed over a period of 
three decades lie within the experimental uncertainties. 

The coupled-cluster calculations, with the exception 
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TABLE XI: Ground and excited (nsnp 3 Po) state scalar polarizabilities Qo (in a.u.) of group II atoms and divalent ions. 
Uncertainties in the last digits are given in parentheses. References are given in square brackets. The abbreviations con- 
form to those used in Table IIVI and [V] Hybrid-RCI+MBPT include experimental data for some transitions. "RCCSDT, 
6 MBPT, c Hybrid-RCI+MBPT data for the alkaline-earth ground states from [57| are listed as recommended "sum-rule" data, 
d RCI+MBPT. 



1 > . . 

J3t? 


IVlg 




C,. 

or 


j_>a 




ivictnoQ 


2s 2 


3s 2 


4s 2 


5s 2 


6s 2 


7s 2 




37 755 [1211 












Th. ECG 


37.73(5) a [226] 


71. 7 b [227J 


157 6 [227] 








Th. 


37.807 [228] 


70.90 [229J 


171.7 [230J 








Th. CI 


37.29 [104] 


70.74 [104] 


156.0 [104] 








Th. CICP 


37.69 [105] 


71.35 [105] 


159.4 [105] 


201.2 [105] 






Th. CICP 






158.00 [35] 


198.85 [35] 


273.9 [35] 


248.56 [35] 


Th. RCCSDT 






152 [133] 


190 [133] 


275.5 [134] 




Th. RCCSDT 


37.76 [57J 


71.33 [57] 


159.0 [5T] 


202.0 [57] 


272.1 [57] 




Th. RCI+MBPT 






169(17) [68] 


186(15) [69] 


268(22) [69J 




Expt. EH 




74.9(2.7) [231] 


157.1(1.3) [57] 


197.2(2) [57] 


273.5(2.0) [57] 




Sum-rule c 


Z D Ap JQ 


QqQt, 3 p° 
o&op rg 


4e4n 3 PS 


5e5n 3 R? 


fiefin 3 R? 






QQ AO TOO] 


1U1.0 IZoz, ZooJ 


one o roo/)l 








in. oiLy-r 














Th HvhriH "RPT-J-MRPT 

_L 11. 1 1 V Ui lu l\v_y 1 1V1 1J JT _L 








^0 t .U 1 1^1] 






± 11. JA 1 | 1V1 D 1 _L 


A1+ 


Si 2 + 


Zn 


Cd 


Hg 


Yb 




3s 2 


3s 2 


4s 2 


5s 2 


5d 10 6s 2 


4d 14 6s 2 




24.2 6 [227] 












Th. 


24.14(12) [235] 


11.688 [88] 


38.12 [113] 


44.63 [113] 


31.32 [113] 




Th. CICP 


24.12 CI [229] 


11.75 CI [229] 
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67.69 [113] 


75.29 [113] 


55.32 [113] 
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54.6 [145] 


315.9 [143] 


Th. RCI+MBPT 



252(25) [240] Th. RCI+MBPT 
266(15) [13] Th. RCI+MBPT 
302(14) [143] Th. RCI+MBPT 



of the RLCCSD one [56[, tended to give polarizabilities 
which were 1 — 2% larger than experiment until the most 
recen t RCCSDT calculation of Thakkar and Lupinetti 
[HI which gave 162.9(6) a.u.. The earlier CCSDT calcu- 
lations tend to overesti mate the polarizabilit y most likely 
due to basis set issues [nl Gil HH, Hi]. The same 
problem could also be leadi ng to the overestimati on of 
the polarizability by the CI and CEPA-PNO pill 
calculations. 



By way of contrast, the RLCCSD calculation (5(| gave 
a polarizability of 163.0 a.u. which is in agreement with 
experiment. We have discussed the differences of the RL- 
CCSD approach from the other coupled-cluster calcula- 
tion in Section TlVFI An important feature here is that 
this calculation uses a i?-spline basis which is effectively 
complete [13, HfJ. As we ha ve discussed on the example of 
He CCSD calculation [147| , polarizability coupled-cluster 
results vary significantly with the choice of the basis set 
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TABLE XII: Ground state polarizabilities ceo (in a.u.) of other systems. Uncertainties in the last digits are given in parentheses. 
References are given in square brackets. The average over magnetic projections is given for atoms which do not have a spherically 
symmetric ground state. "Third-order MBPT, 6 RRPA, C RCCSD, d RLCCSD, e CICP, / light deflection, 9 RCI+MBPT, h RESIS 
reanalysis using theoretical estimates of higher order polarization corrections, 'spectral analysis, -'CI, fc CCSDT, 'f-sum rule, 
m RMBPT 
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57.74 fe [36] 






49.2 9 [251] 






Th. 


59.5 j [252] 












Th. 




49.9 [253] 


61.9 [253] 


51.6 [253] 






Th. RCCSDT 


46.2(20) [2541 




68.7(8.1) [2551 


51.3(5.4) [256] 






Expt. EH 


Si 


Sn 


Pb 


Ir 
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3s 2 3p 2 


5s 2 5p 2 


6s 2 6p 2 










37.0 [252] 












Th. CI 


37.17 [36] 












Th. CCSD 


37.3 [257] 


52.9 [257] 


47.3 [2571 








Th. RCCSDT 




42.4(11.0) [257] 


47.1(7.0) [257] 


54.0(6.7) [258] 


137.0(9.4) / [76] 




Expt. EH 



if it is not sufficiently saturated. In summary, large (ef- 
fectively complete) basis sets are needed for precision po- 
larizability calculations by a coupled-cluster method. 

The relativistic corre ction to the dipole polarizability 
is about —1.0 a.u. |l74j . The three non-relativistic CICP 
calculations all lie within 0.5% of the experimental polar- 
izability. As mentioned earlier, these calculations implic- 
itly include relativistic effects to some extent by tuning 
the core polarization potential to the experimental bind- 
ing energies. The RLCCSD calculation uses a relativistic 
Hamiltonian and intrinsically includes relativistic correc- 
tions. The recommended value 162.6(3) is based on the 
RLCCSD calculation with resonant 3s — 2>pj transition 
matrix elements replaced by their experimental values. 

The Na polarizability of 162.7(8) [72[ obtained by in- 
tcrfcrometry experiment served as the reference polariz- 
ability in the determination of the K and Rb polarizabil- 
ities by the interferometry ratio approach [73] . Table IVl 
shows excellent agreement of these values with the hybrid 
RLCCSD f-sum polarizabilities of [56| . 



3. Mg+ 

We use Mg + to illustrate the RESIS experimental ap- 
proach owing to recent advances in that area. Both the 
potential and the problems of determining the polariz- 
abilities of ions using spectral analysis are evident by 
contrasting the different values listed for Mg + and Ba + . 
The original analysis of the RESIS data f or M g + reported 
a dipole polarizability of 35.00(5) a.u. 198]. However, 
the contributions from the Cj and C$ terms of Eq. (|3"6")l 
can possibly corrupt the value of «o if they are significant 
as described in Section [ill II A more detailed analysis of 
the polarization plot which explicitly included the C7 and 
Cs terms was subsequently performed in Ref. [84j. This 
polarization plot is shown in Figure [2] The data points 
including the explicit subtraction of the C7 and Cs terms 
show a higher degree of linearity. The revised analysis re- 
sulted in ao — 35.05(3) a.u.. This is only 0.15% larger 
than the original value and lies within the original error 
limits. 

The treatment of non-adiabatic corrections is a much 
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TABLE XIII: Excited state scalar ao and tensor a 2 polariz- 
abilities (in a.u.) of selected systems. Uncertainties in the 
last digits are given in parentheses. References are given in 
square brackets. a CICP, 6 RCI+MBPT, C CCSDT, d RCCSDT, 
e one-electron model potential. 
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242. 4 a [234] 






a ■> 
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278(15) 6 [240J 
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24.26(84) 
23.35(52) 
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more serious issue for the Ba + ground state. Table IVl 
show s that subsequent analysis of the RESIS data [l97l 
Il99l ] do not lie within their mutual uncertainties. The 
most recent analy sis o f RESIS data gave a polarizability 
of 123.88(5) a.u. |l99l ]. This analysis explicitly included 
non-adiabatic effects from the low-lying bd excitation. 
However, non-adiabatic effects from the 6s-6p excitation 
are also significant and need to be included for a RESIS 
polarizability to be regarded as definitive. 

The influence of the non-adiabatic effects in the CV and 
Cs terms of Eq. (|36[) can be minimized by taking mea- 
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FIG. 2: The polarization plot of the fine-structure intervals 
of Mg for the n = 17 Rydberg levels. The A_B c i intervals are 
corrected for relativistic, second-order and Stark shifts. The 
AE c2 intervals account for (r -7 ) and (r -8 ) shifts. The linear 
regression for the AE c2 plot did not include the last point. 



surements at high values of L, e.g. L > 8. Unfortunately, 
as the non-adiabatic corrections diminish with increasing 
L, the states with very high L are more sensitive to Stark 
shifts due to stray electric fields. As the energy splitting 
of the Rydberg states gets smaller at higher L, the po- 
larizabilities of the (n, L) levels then get larger due to 
the very small (n, L — n, L ± 1) energy differences. To 
a certain extent one has to choos e th e (n,L) states to 
navigate between the low-L Scylla [286[ of non-adiabatic 
corrections and the high-L Charybdis 286] of Stark shifts 



4- Cesium 

The Cs atom has been studied extensively o wing to 
the parity- violation experiments on this system [287| . A 
comprehensive set of Cs scalar and tensor polarizabilities 
for the 7s - 12s, 7p 1/2 - 10pi/ 2 , 7p^/ 2 - 10p 3 /Q, 5«f a 



lOd 



12, 



and 5d 



5/2 



lQd fi 



states 
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177, 20' 



I2ia4225| taken from Ref. [n3 is given in Table E 

The polarizabilities listed in Table |X] are in 10 3 a.u. 
since the values range in size from 300 a.u. to 7 x 10 6 
a.u. The results of Ref. jl77| are obtained from sum- 
over-state calculation using the RLCCSDT matrix ele- 
ments and experimental energies for a large number of 
states. The remaining contributions from highly-excited 
states were evaluated as well. In a few cases, some of 
the RLCCSDT matrix elements have been r eplac ed with 
matrix elements extracted from experiment [1771 ] . Incor- 
porating such highly-excited states as 12s required the 
use of a very large R — 220 a.u. spherical cavity and 
large B-spline basis sets. Extensive tests of numerical 
stability of the calculations in such a large cavity have 
been conducted to verify the accuracy of a finite basis 
set representation. All matrix elements used to evaluate 
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dominant polarizability contributions were critically eval- 
uated for their accuracy based of the size and type of the 
dominant correlation corrections and semi-empirical esti- 
mates of the omitted correlation terms. Such uncertainty 
evaluation is discussed in more detail in Section [VII 

Coulomb approximation (CA) values |218| were also 
computed with a sum-over-states approach. One inter- 
esting feature of Table [X] is the reasonable level of agree- 
ment between the CA and RLCCSDT values for many of 
the polarizabilities. The CA results are computed with 
wave functions which are tuned to experimental energies. 
The radial matrix elements that arise in the sum-over- 
states calculation are dominated by the form of the wave 
function at large distances. Tuning the wave functions to 
have the correct energy goes a long way to ensuring that 
the long-range part of the wave function has the correct 
shape. 

A number of the experimental values in Table IXl were 
obtained from Stark shift experiments. In many cases, 
the excited state polarizabilities are much larger than 
the Cs ground state polarizability, so uncertainties in the 
ground state have minimal impact on the overall uncer- 
tainty. The agreement between the experimental and RL- 
CCSDT polarizabilities is excellent for the ns states, in 
most cases the difference between them is less than 1%. 
The situation is not so clear-cut for the nd states. Differ- 
ences between theory and experiment are large in some 
cases, but so are the uncertainties of many of the ex- 
perimental values. However, the RLCCSDT results were 
found in good agreement with more recent ex perim ents 
[2221 [Ha SH- The RLCCSDT calculation [177| pro- 
vided critically evaluated recommended values for a large 
number of Cs polarizabilities for which accurate experi- 
mental data are not available. 



C. Two electron atoms and ions, ns 2 1 S and 
nsnp 3 Pq states 

Table IXII gives the polarizabilities for a number of di- 
valent species including the alkaline-ear th a t oms from 
Refs. r 




„ 28, \M 153|, p , | 6 4 16 8 1 16 ft M 10 4 1051 . 11131 
Il2ll . fl33Tll3l . Il7lll43lll^^ [236rl240[ EH Exf . 

The beryllium atom serves as a theoretical benchmark 
since a very accurate value has been obtained with aba- 
sis of exponentially correlated Gaussians (ECG) |l2l| . 
The CICP [l0| and RCI+MBPT polarizabilities [HJlie 
within 0.2% of the ECG basis polarizability. 

The sub-1% agreement between the highest quality 
theory and experiment that occurred for the alkali atoms 
is not observed for the alkaline-earth atoms owing to 
their more complicated atomic structure and resulting 
mixing of configurations. As we have described in Sec- 
tion IIV G[ perturbative methods do not work well for 
strong valence- valence correlations. The hybrid val- 
ues for Ca and Sr based on the RCI+MBPT calcula- 
tions with the matrix elements for the resonance tran- 
sitions replaced by values derived from experiments are 



TABLE XIV: Static polarizability differences (in a.u.) de- 
rived from selected Stark shift measurements. Uncertainties 
in the last digits are given in parentheses. References are 
given in square brackets. "RMBPT, ^RLCCSDT, C CICP, 
d RCI+MBPT. 



Atom 


State 
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Theory 
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6s — 7s 


5837(6) f93] 


5834° [288] 






5709(19) 


[289] 






6s-8s 
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[216] 


37820(290) 6 [216] 
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32. f (4.0) 


[290] 


37.6 C [232J 
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16.3 C [232J 
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(m = 0) 


90.4(13.5) 


[2921 


107.5 C [234J 






98.97(33) 


[2931 




Ba 


6s 2 - 6s6p 1 Pi 










(m = 0) 


-229.32(48) 


[294] 


-247 d [53] 


Yb 


6s 2 - 6s6p 3 P° 






160(60) d [240] 




(m = 0) 


123.85(38) 


[278] 


110(18) d [240J 


Hg 


6s 2 - 6s6p 3 P° 


26.68(48) 


[295] 


26.95 d [145] 


Ga 


4s 2 4p 3/2 - 4s 2 5s 


788(40) 


[183] 




Tl 


6s 2 6pi/ 2 - 6s 2 7s 


-900(48) 


[296] 


-830 6 [297] 






-829.7(3.1) 


[298] 




Tl 


6P1/2- 7pi/a 


-4967(249) 


[296] 


-4866 b [297] 


Yb+ 


6s-5d 2 D 3/2 


-41.8(8.5) 


[279] 





respectively 1.1% and 2.5% smaller than the ab initio 
RCI+MBPT estimates [57} • With the exception of Be, 
our recommended values for alkaline-earth polarizabil- 
ities are those obtained from the hybrid RCI+MBPT 
method. We note very good agreement of the RCCSDT 
calculations of Ref. [35[ for the ground state polarizabil- 
ities of Ca, Sr, and Ba with the recommended values 
in all three cases. One of the problems of the hybrid 
approach is the paucity of high-precision experimental 
data for divalent atoms. Strontium is the only atom 
where the polarizability has been quoted with a preci- 
sion approaching 0.1% [57j • This is due to the availabil- 
ity of a high precision estimate of the resonant oscillator 
strength obtained by Yasuda and Katori using photo- 
Howe ver, an alternate 
gave a lifetime 0.8% 



association spectroscopy [301 1 



photo-association experiment [302| 
smaller than the Yasuda and Katori value, so it may be 
over-optimistic to assign an uncertainty of 0.1% to the 
strontium polarizability. Currently the best estimate of 
the 5s5p 3 Pq excited state polarizability of Sr is accurate 
to 0.8% despite the use of the experimental da ta. T he Sr 
polarizabilities are discussed in detail in Ref. jl4l| . 

The 27 A1 + ion is included in Table IXII since it is be- 
ing used in the d evelo pment of a single ion optical fre- 
quency standard 



303] 



The most reliable calculation of 
the ground state polarizability ao is probably given by 
the CICP calculation. The only experimental value is of 
low precision (3%) and was obtained by summing exper- 
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imental oscillator strengths. A CICP calculation of the 
isoelectronic Si 2+ system gave a polarizability that was 
within 0.2% of the value from a RESIS experiment. 

The scatter amongst the different calculations of yt- 
terbium underlines the difficulties of performing calcula- 
tions in this system. The source of the problem lies in 
the weakly bound 4/ 14 core. There are 20% differences 
between two of the RCI+MBPT calculations that are 
discussed in recent work by Dzuba and Derevianko [l43j | 
and are attributed to the inconsistent use of experimen- 
tal matrix element for the principal transition in [f3j . Yb 
is of particular interest for many applications, including 
ultracold atoms, optical frequency standards, and parity 
violation experiments. 

There is a significant discrepancy for Cd between the 
refractive index value of 49.65(f .49) a.u. [64| and the cal- 
culated value of 44.63 a.u. from the RCICP calculation 
[113| . For a number of reasons, including the measured 
values of the oscillator strengths for the 5s 2 1 S f — 5s5p 
1 P° transitions, it has been suggested that t he e xperi- 
mental polarizability might be overestimated |304| . 

The polarizabilities of other excited states, tensor po- 
larizabilitics, and Stark shifts in divalent systems are dis- 
cussed in the next subsection. 



D. Other data 

Ground state polarizabilities for the other selecte d sys - 
tems from Refe. B M , M, Hi EM [HI [24lT[245l l247i - 
25il . l25fl . [2521 . [252^581 [305^ are given in Table EH In 
this review, we list data for selected systems with a mono- 
valent ns ground state: Cu, Ag, Au, Zn+, Hg+, and 
Yb + , and Al 2+ ; ions for which recent RESIS experiments 
have been performed: Si 3+ and Kr 6+ ; neutral atoms with 
three and four valence electrons: Al, Ga, In, Tl, Si, Sn, 
Pd. and Ir; and U. The reader is referred to a recent re- 
view [2l| for atomic ground state polarizabilities of other 
systems not listed herein. 

One notable discrepancy between theory and experi- 
ment occurs for the Al ground state where the best calcu- 
lations exceed the experi men t value from an EH balance 
experiment by 25% [36|, |254| . The most precise experi- 
mental value in Table |XTI1 is the RESIS value for Si 3+ . 
The final value, 7.433(25) a.u. comes fro m a reanaly- 
sis of the raw experimental data [248l l306l | that includes 
estimates of r~ 7 and r~ 8 polarization corrections from 
RLCCSDT and CICP calculations 



Th e ag reement 

between the RLCCSD polarizability of 7.419 [Hj and the 
latest RESIS reanalysis is at the 0.2% level. 

Table lXllTl shows a number of measurements and calcu- 
lations of the tensor polarizability of non-alkali systems 
including Ca Jsr, Ba, Zm Cd, Hg, Tl, Yb, Y b+ from 

Refs m. IBl. jiiiL pa \ml [23I I23I I240I I25I I253I [259! 

[2591 . [259M2661 I268L I270U273L I275M279L l293| These sys- 
tems are the ones under consideration as frequency stan- 
dards or are being used in atomic parity violation experi- 
ments. Measurements for some states have been omitted 



from the Table, and some older or les s pre cise results on 
Sc, Y , La and Lu [307|, Cd 206], Ba Hi HM Hoi , Hg 
[206| . Yb [2761 ], Sm and Eu [273] have also been omitted. 

One feature of Table IXIIII is the relatively small num- 
ber of modern calculations performed. For example, the 
best calculated polarizabilities for the 4s4p 1 > 3 p° states 
of Ca are the non-relativistic CICP calculations. An- 
other feature is the relatively large uncertainties in many 
of the experimental values. There are only five tensor po- 
larizabilities with uncertainties less than 2%. The most 
precisely measured a2 of —43.04(40) a.u. occurs for the 
Ba 6s6p l P? state. The RCI+MBPT value of -51 a.u. 
is incompatible with experiment. 

The static polarizability differences (in a.u.) for se- 
lected transitions in Cs, Mg, Ca, Ba, Yb, Hg, Ga . Tl , 
and Yb + derived from Stark shif t mea surements [93lll83l 

[M IM IM IMS IIM 

are compared with 



theor etical calculations [5a Il45l 12161 [2341 . [28l [2971 ] in 
Table IXIVl Total polarizability differences are given for 
the cases where m values are listed, otherwise scalar po- 
larizability differences are listed. 

There have been sub-1% experiments on four systems, 
Cs, Ba, Yb, and Hg. The ability of RCI+MBPT calcu- 
lations to reproduce experiment for the divalent systems 
is mixed. The agreement for the Hg 6s 2 1 S' - 6s6p 3 P° 
is excellent, but 10% discrepancies exist for the Ba 6s 2 
1 S - 6s6p 1 P? and Yb 6s 2 1 S - 6s6p 3 P? transitions. 
However, the RCI+MBPT calculations for Ba [53[ and 
Yb [HI were among the first RCI+MBPT calculations 
reported. 



VI. EVALUATING UNCERTAINTIES OF 
THEORETICAL VALUES 

A. Sources of theoretical uncertainty 

As illustrated by the tables in the previous section, 
benchmark comparisons of theory and experiment carry 
more value when the theoretical results are accompanied 
by uncertainty evaluations. Uncertainty bounds are par- 
ticularly important for the recommended values obtained 
by either high-precision theory methods or by combina- 
tion of theory values with experimental data. The analy- 
sis of the theoretical uncertainties has been stimulated by 
the applications that require an error bound to be placed 
on the recommended values. Such applications include 
parity violation, development of the next-generation fre- 
quency standards, ultra-cold atom studies, etc. Analysis 
of certain experiments requires input of some data that 
cannot be easily measured and have to be obtained from 
theory. In those cases, the uncertainties of the theoret- 
ical data have to be included in the uncertainty of the 
final experimental value. Evaluations of the theoretical 
uncertainties are still few and cannot be carried out for 
all of the methods and in all cases. Here, we discuss how 
some theoretical uncertainties may be evaluated. 

There are two distinct sources of theoretical uncertain- 
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ties. First, there is an uncertainty associated with the 
numerical constraints upon the calculations. Many of 
the methods that we discussed in this review are com- 
putationally very intensive and restrictions are imposed 
so that the calculations can be performed within a rea- 
sonable time. Most common numerical uncertainties are 
associated with the choice of the basis sets, configuration 
space, radial grid, termination of the iterative procedures 
after achieving the specified convergence tolerance, etc. 
Generally, it is possible to at least estimate uncertainties 
caused by numerical issues by varying the appropriate 
parameters and recording the changes in the results. In 
many cases, it is possible to simply continue to change 
parameters until the change in the resulting values is suf- 
ficiently small or negligible. 

For example, it is relatively easy to test the conver- 
gence of B-spline basis sets. The dimensionality of the ra- 
dial basis in a RLCCSD calculation for each partial wave 
(e.g, ns, npi/2, np 3 / 2 > ■■■ states) is steadily increased. 
The final values of a property like the sodium ground 
state polarizability do not change, within the quoted dig- 
its whether the B- spline basis has a dimension of 40, 
50, or 70 orbitals [56[. Using only 20 orbitals, however, 
will lead to change in the final value that is not negli- 
gible. Also, truncating partial wave expansion at I = 3 
will measurably affect the final result, while including all 
partial waves up to I = 6 is sufficiently complete in this 
case. Generally, such tests do not have to be carried out 
at the level of the most accurate calculation possible and 
it is sometimes sufficient to study the lowest-order re- 
sults or low-order MBPT values. In some cases, it may 
become necessary to completely repeat the entire calcu- 
lation. However, such numerical problems may be stud- 
ied by well understood conventional methods. In most 
cases, numerical errors of the theoretical values can be 
made small enough not to affect any of the significant 
figures that are quoted, or can be evaluated and quoted 
as uncertainty in the last digit. 

Investigations using the Hylleraas method typically 
perform a series of calculations of increasing dimension 
while keeping the non-linear parameters the same. The 
convergence of the data against a value of the total poly- 
nomial power is studied. The total polynomial power for 
a correlated wave function such as Eq. (|43p would be 



^ = h + h + h + hi + ji3 + h3- 



(63) 



Most expectation values in a Hylleraas calculation con- 
verge as ~ l/il p . This result is exploited to give uncer- 
tainties in energies, transition matrix e lements, polariz- 
abilities, and other quantities [ll7l Ill8j . 

The theoretical uncertainties of the second type are 
much harder to evaluate. These are the uncertainties as- 
sociated with the particular theoretical methodology, for 
example, the uncertainty associated with stopping a per- 
turbation theory treatment at third order. Ideally, the 
total uncertainty of the theoretical value should give an 
estimate of how far any value is from the actual (un- 
known) exact result. Evaluation of the complete theoret- 



ical uncertainty is non-trivial since it essentially involves 
the evaluation of a quantity that is not known beforehand 
and cannot be determined by the theoretical methodol- 
ogy adopted. 



B. Sources of uncertainties in the sum-over-states 
polarizability calculations 

It is particularly problematic to evaluate full theoreti- 
cal uncertainties for the semi-empirical theoretical meth- 
ods. In this case, there may be no basis to make assump- 
tions regarding the missing theory. It may be possible to 
infer some information based on the agreement of CICP 
calculations with quality experiments for similar states 
in other members of the same iso-electronic series. For 
example, the CICP ground state polarizability for Al + of 
24.14 a.u. has been assessed at ±0.5% [235j ] based on the 
0.3% agreement between a CICP calculati on of the Si 2+ 
polarizability and a RESIS experiment [HI, |239| . The as- 
sessment of uncertainties, for states that lack validating 
information, as in the case of th e 3 Pg state of Al + has a 
larger speculative element 235 1. 

Several strategies exist for uncertainty evaluation 
for the ab initio MBPT, correlation potential, and 
all-order linearized coupled-cluster (RLCCSDT) ap- 
proaches. These strategies are illustrated using RL- 
CCSDT method which utilizes the sum over states al- 
gorithm. For brevity, we refer to RLCCSDT calculation 
as "all-order" in the text below. 

We use the example discussed in Section lH A 21 i.e. the 
polarizability of the 5pi/2 state. Table Hill lists a detailed 
breakdown of the contributions to this value. There are 
three separate contributions: the main part (5s — lis 



and 4^3/2 — 9^3/2), remainder (all other valence terms), 
and core contribution. The uncertainty in each term of 
the main part has to be determined. The energy lev- 
els of low-lying states are generally well known. There- 
fore, the determination of the uncertainty here reduces 
to the evaluation of the uncertainty in the corresponding 
electric-dipole matrix elements. The relative uncertainty 
in the polarizability contribution is twice the relative un- 
certainty in the matrix element (see Eq. (|17[) ). 

The uncertainty of the remainder (higher n contribu- 
tions) as well as the uncertainty of the ionic core have to 
be determined separately. The uncertainty in the RPA 
value of the core is estimated from comparison of the 
RPA values for noble gases with experiment and preci- 
sion coupled-cluster calculations (see Table IIVI and the 
corresponding discussion). The evaluation of the uncer- 
tainty of the remaining highly-excited contribution has 
been discussed in gr eat detail in recent work on the Sr + 
polarizabilities [193|. 

In most cases, all of the uncertainties are added in 
quadrature to obtain the final uncertainty of the polar- 
izability value. 
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1. Determination of the uncertainties in El matrix 
elements 

Ultimately, the theoretical uncertainty estimates in the 
polarizability need uncertainties in the El matrix ele- 
ments such as those listed in Table IIIII The starting 
point of relativistic MBPT or all-order RLCCSD calcu- 
lations for monovalent systems is a DHF calculation. We 
refer to the DHF value as the lowest order. Essentially 
all corrections to that value come from Coulomb corre- 
lations. Breit interaction corrections to the El matrix 
elements are g enerally insignificant at the present level of 
accuracy [3091 ] , and the relativistic corrections are intrin- 
sically included due to use of a relativistic Hamiltonian. 
Therefore, an uncertainty evaluation requires an estima- 
tion of the missing part of the correlation correction. The 
strategies to do so include: 

• approximate evaluation of the size of the correlation 
correction; 

• evaluation of the size of the higher-order correc- 
tions; 

• study of the order-by-order convergence of pertur- 
bation theory; 

• study of the breakdown of the various all-order con- 
tributions and identification of the most important 
terms; 

• semi-empirical determination of dominant missing 
contributions. 

The first three strategies are aimed at providing rough 
estimate of the matrix element uncertainty. Separate 
third-order RMBPT and all-order calculations have to 
be carried out to evaluate the accuracy of the all-order 
values since the extraction of third-order matrix elements 
from the all-order values is impractical. 

The application of the first three strategies are illus- 
trated in Table |XV] where Rb El matrix elements are 
listed [310( | . Three values are given for each matrix ele- 
ment: lowest-order DHF value, the third-order RMBPT 
value, and all-order values obtained from an RLCCSD 
calculation. Third-order values include the second-order, 
third-order, and RPA corrections iterated to all-orders 
(see [ll| for detailed description of the third-order MBPT 
calculations). The size of the correlation correction is 
estimated as the relative difference between the lowest- 
order and the all-order values. It is given as a percentage 
change in the rows labelled "correlation" . The size of the 
fourth and higher-order corrections is estimated as a per- 
centage difference between the third-order and all-order 
values and listed in the rows labelled "higher orders" . 

Study of the "correlation" and "higher orders" rows 
gives some insight into the accuracy of the final all-order 
values. First, it is noted that the corrections vary sig- 
nificantly among the different transitions. Very rough 
estimate of the uncertainty can be obtained by assuming 



TABLE XV: Rb electric-dipole matr ix el ements (in a.u.) cal- 
culated in different approximations [310(. The rows labelled 
"correlation" list an estimate of the correlation contribution, 
determined as the relative difference between the lowest-order 
and the all-order values. The rows labelled "higher orders" list 
an estimate of the 4th and higher-order contributions, deter- 
mined as the relative difference between the third-order and 
the all-order values. Absolute values are listed. The negative 
sign in front of the lowest-order 6^3/2 — Qpi/i value indicates 
that the lowest-order gives incorrect sign for this matrix ele- 
ment. 





5s - 5pi /2 


5s - 6pi /2 


6s - 5pi /2 


Lowest order 


4.819 


0.383 


4.256 


Third order 


4.181 


0.363 


4.189 


All order 


4.221 


0.333 


4.119 


Correlation 


14% 


15% 


3.3% 


Higher orders 


0.9% 


9% 


1.7% 




8s - 8pi /2 


4d 3/2 - 5pi /2 


6d 3/2 - 6pi /2 


Lowest order 


26.817 


9.046 


-0.047 


Third order 


25.587 


8.092 


2.184 


All order 


25.831 


7.847 


2.974 


Correlation 


3.8% 


15% 


100% 


Higher orders 


0.9% 


3% 


27% 



that higher-order corrections incorporated into RLCCSD 
are smaller than the higher orders that are omitted by 
RLCCSD. Thus, the difference between third and all- 
orders is taken as the uncertainty. In most cases, this 
procedure will significantly overestimate the uncertainty 
since Table IXVI shows that contributions from all higher 
orders are lower than the second and third-order in all 
cases except the small 5s — Qpi/2 matrix element. How- 
ever, this procedure clearly indicates that while 5s—5p±/ 2 
matrix element is probably accurate to better than 1%, 
SD all-order 6d 3 / 2 — 6p 3 / 2 matrix element may only be 
accurate to about 25%. 

The last two strategies should be employed if more 
accurate uncertainty evaluations are required. This can 
only be done for certain cases within the framework of 
the RLCCSDT method and requires substantial addi- 
tional calculations and the careful analysis of all avail- 
able data. First, the breakdown of the all-order terms 
have to be studied. Triple excitations need to be added 
at least partially. If certain types of the contributions 
(associated with so-called Brueckner orbital terms) are 
dominant, they may be estimated by the semi-empirical 
scaling described, for example, in Ref. [13] ■ This proce- 
dure involves rescaling single-excitation coefficients p mv 
(see Section HVF[) using experimental energies, and re- 
running the entire matrix element calculation with the 
modified coefficients. Obviously, this method is only ex- 
pected to produce more accurate values if correlation 
correction is dominated by the terms containing single 
valence excitation coefficients. However, this is true in 
many cases. Non-linear terms may also be evaluated. 
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The most extensive uncertainty study of this type has 
recently been carried out for the atom ic quadrupole mo- 
ments of Ca + , Sr + , and Ba + in [31l| . 

Detailed studies of the uncertainties of the electric- 
dipole m atrix elem ents are described, for example, in 
Refs. [30L Il93l l216| . A brief description is given he re fo r 



the case of the 4d 5 / 2 — 5p 3 / 2 matrix element in Sr + [1931 ]. 
This transition is important in the evaluation of the Sr + 
BBR shift. Correlation corrections change the matrix el- 
ement by about 20%. Study of the correction breakdown 
indicates that the correlation is dominated by a single 
term that contains single valence excitations. Therefore, 
we carry out additional ab initio calculations that par- 
tially include triple excitations, and also perform scaled 
RLCCSD and RLCCSDT calculations. The results of 
these four calculations are listed below. All data are in 
atomic units. The first line corresponds to the "all-order" 
lines in Table [XV] 



RLCCSD 
RLCCSDT 
RLCCSD scaled 
RLCCSDT scaled 
Final 



4.150 
4.198 
4.187 
4.173 
4.187(14) 



Note that scaled values are much closer together than 
the SD and SDT ab initio values. The final value was 
taken to be RLCCSD scaled 4.187(14) result (see, for 
example, [37[ and references therein for the discussion of 
this choice). The uncertainty of 0.014 is determined as 
the maximum difference between the scaled SD values 
and the ab initio SDT and scaled SDT values. 



VII. APPLICATIONS 

A. Parity non-conservation 

The goals of the parity nonconservation (PNC) studies 
in heavy atoms are to search for new physics beyond the 
standard model of the electroweak interaction by precise 
evaluation of the weak charge Q Wl and to probe parity 
violation in the nucleus by evaluating the nuclear anapole 
moment. The study of PNC in the cesium 6s — 7s tr ansi- 
tion involving both high-precision measurement [287| and 
several high-precision calculations provided an atomic- 
physics tes t of the standard model of the electroweak in- 
teractions 312]. Moreover, accurate determination of the 
uncertainty in theoretical values was necessary, leading to 
detailed studies of parity-conserving quantities in Cs in- 
cluding the polarizab ilities of the 6s, 6pj, and 7s states 
(see [3jJ, |37J, I313L l314j ] and references therein) . The anal- 
ysis of the Cs experiment was instrumental in developing 
methods to evaluate the uncertainties of the theoretical 
data [315j ]. 

In the Cs experiment [287j . the PNC amplitude was 
measured relative to Stark-induced tensor transition po- 



larizability j3s (some works refer to this quantity as the 
vector transition polarizability). The DC electric field 
mixes states of opposite parity allowing electric-dipolc 
transitions between ns states. The Stark- induced ampli- 
tude is expressed via the Stark-induced scalar and tensor 
transition polarizabilities as and (3 s- In the case of the 
Cs 7s — 6s transition, they are calculated as sum-over- 
states using the expressions [3151 ] 



«s = j ^2 ( 7s \\ D \\ n Pi/2)( n Pi/2\\D\\6s 



lj2( 7s \\D\\np 3/2 ){np 3 / 2 \\D\\6s) 



1 1 

+ 



Et s — E np3/2 Eq s — E. 



nps/2 



^ = ^E< 7S ll^ll^l/2>(^l/2ll^l|6s) 



E-js E npi/2 Eq s E npi/2 J 
\j2(7s\\D\\np 3/2 )(np 3/2 \\D\\6s 
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E-js E np3/2 Eq s E np3/2 



(64) 



These quantities have been extensi vely studied d ue to 
their importance in PNC research, [287l 1313143161 . It is 
more complicated to calculate (3s accurately, in compar- 
ison to as, owing to severe cancellations between dif- 
ferent terms contributing to (3 s- The ra tio o f as and 
fis has been measured to high precision [316| . At the 
present time , Cs experiment is consistent with the stan- 
dard model [312| . 

Ho wever , the precise measurement of PNC amplitudes 
in Cs 287] also led to an experimental value of the small 
contribution from the nuclear-spin dependent PNC ac- 
curate to 14%. The constraints on weak nucleon-nucleon 
coupling constants derived from this experiment and cal- 
culations in Cs were found to be significantly inconsistent 
with constraints fro m deep ine lastic scattering and other 
nuclear experiments [317H319| . At the present time, this 
discrepancy remains unexplained. 

More PNC experiments in other atomic systems, such 
as Ra + , Yb, and Fr are currently in progress. Experi- 
ments in Pb, Bi, and Tl have been conducted but theoret- 
ical calculations of comparable accuracy are not available 
to permit precise comparison of experiments with the 
standard model. Comparison of theoretical and e xperi - 
mental values of as and (3s for Tl are given in Ref. [2971 ] . 
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B. Ultracold atoms in optical lattices and quantum 
computation 

Quantum computation is a field of research which is 
aimed at using the quantum nature of matter to pro- 
duce fundamentally new methods of computation. There 
are various approaches to the experimental realization of 
the quantum computation. In the quantum computation 
scheme relevant to this review, the qubits are realized 
as internal states of neutral atoms trapped in optical lat- 
tices or microtraps. This approach to quantum computa- 
tion has many advantages, such as the long decoherence 
times of the internal states of the atoms, flexibility in con- 
trolling atomic interactions, scalability, possible massive 
parallelism, and well-developed experimental techniques. 

Trapping an atom or a group of atoms in an optical 
lattice raises the possibility that the laser field used to 
create the lattice might shift the energy levels of the lower 
and upper states by different amounts. This can result 
in a wavelength (and intensity) dependent shift of the 
clock transition. This issue was hrst raised for the atomic 
clocks based on neutral atoms trapped in optical lattices. 

A solution to this problem was proposed by Katori et 
al. [HI who suggested that the laser can be tuned to a 
magic wavelength A mag i c , where lattice potentials of equal 
depth are produced for the two electronic states of the 
clock transition. At such wavelength, the AC polarizabil- 
ities of the two relevant states satisfy the condition 

Crupper (A) — Ot lower (A). (65) 

Ah initio calculations of the dynamic polarizability are 
valuable in making an initial estimate of the magic wave- 
length prior to construction of the optical lattice. How- 
ever, it is possible to make very precise determinations 
of the magic wavelength once the lattice has been con- 
structed and atoms have been trapped since the exper- 
imental design is that of a null experiment. The exper- 
imental magic wavelength can be used as a constraint 
upon the dynamic polarizability and used to refine the 
polarizability calculation. 

Examples of a magic wavelength calculation are de- 
picted in Fig. |3l where polarizabilities of the Li 2s and 
3pi/ 2 states obtained using RLCCSDT method are plot- 
ted. The magic wavelengths are located at the crossing 
points of the two curves. The ground state polarizability 
is nearly flat in this wavelength region, while the 3pi/2 
polarizability has several resonances noted by the vertical 
lines. 

One of the current goals of the quantum informa- 
tion projects is to design an apparatus capable of in- 
terconnecting "flying" and "stationary" qubits. The 
ability to trap neutral atoms inside high-Q cavities in 
the strong coupling regime is of particular importance 
for such schemes. In a far-detuned optical dipole trap, 
the potential experienced by an atom can be either at- 
tractive or repulsive depending on the sign of the AC 
Stark shift due to the trap light. The excited states 
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FIG. 3: Magic wavelengths for the 2s — 3pi/2 transition in 
Li. Upper state of the resonant transition is marked on top 
of the box. 



may experience an AC Stark shift with an opposite 
sign to the ground state Stark shift which will affect 
the fidelity of the experiments. McKeever et al. |320| 
demonstrated state-insensitive trapping of Cs atoms at 
Amagic = 935 nm while still maintaining strong coupling 
for the 6P3/2 — 6S1/2 transition. 

The magic wavelengths in Na, K, Rb, and Cs atoms 
for which the ns ground state and either of the first two 
npj excited states experience the same optical potential 
for state-insensitive cooling and trapping were evaluated 
in [29[ . This was accomplished by matching the dynamic 
polarizabilities of the atomic ns and npj states using ex- 
tensive relativistic all-order calculations. Uncertainties 
in the dynamic polarizabilities were also evaluated. 

One requirement for the experimental realization of the 
scalable quantum computer is the design of a quantum 
gate with low error rate which will allow for error cor- 
rection. Therefore, it is important to study the various 
decoherence mechanisms and to search for ways to opti- 
mize gate performance. 

The issue of the mismatch of the polarizabilities of the 
ground and excited states has also arisen in schemes to 
perform quantum logical operations where it i s a s ource 
of decoherence. In the Rydberg gate scheme |32l| . the 
qubit is based on two ground hyperfine states of neu- 
tral atoms confined in an optical lattice. A two-qubit 
phase gate may be realized by conditionally exciting two 
atoms to relatively low-lying Rydberg states. The choice 
of this particular scheme results from its potential for fast 
(sub-microsecond) gate operations. Su ch a gate has been 
recently experimentally demonstrated |322| . An atom in 
a Rydberg state will, in general, move in a different opti- 
cal lattice potential than that experienced by the ground 
state. Therefore, the vibrational state of the atom in 
the lattice may change after the gate operation is com- 
pleted, leading to decoherence due to motional heating. 
The optical potential for a given state depends on its AC 
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polarizability, so we can seek to minimize this motional 
heating effect by the choice of a particular Rydberg state 
or of the lattice photon frequency uj. A method for ac- 
complishing this by matching the frequency-dependent 
polarizabilities a(A) of the atomic gr ound state and Ry- 
dberg state is des cribe d in [323l l324j . 

In recent work [325j , a novel approach to quantum in- 
formation processing, in which multiple qubits can be 
encoded and manipulated using electronic and nuclear 
degrees of freedom associated with individual alkaline- 
earth-metal atoms trapped in an optical lattice, was pro- 
posed and analyzed. In this scheme, curves of dynamic 
polarizabilities are needed for alkali and group II atom el- 
ements to locate the wavelengths where one of the species 
can escape or where AC Stark shifts cancel for a specific 
transition. 



C. Atomic clocks 

The current definition of the second in the Interna- 
tional System of Units (SI) is based on the microwave 
transition between the two hyperfine levels of the ground 
state of 133 Cs [343| . The present relative standard uncer- 
tainty of the Cs microwave frequency standard is around 
5 x 10~ 16 . More accurate clocks are needed for a vari- 
ety of applications. Significant recent progress in optical 
spectroscopy and measurement techniques has led to the 
achievement of relative standard uncertainties in opti- 
cal frequency standards that are comparable to the Cs 
microwave benchmark. The frequencies of feasible opti- 
cal clock transitions are five orders of magnitude greater 
than the standard microwave transitions, and so smaller 
relative uncertainties are potentially achievable. A list of 
optical transitions recommended for this purpose has re- 
cently been disseminated b y the International Committee 
for Weights and Measures [344j - 

There are two types of optical atomic clocks under ac- 
tive investigation at the moment. Both types of clocks 
are based on optical frequency transitions with a nar- 
row linewidth. The narrow linewidth mandates that 
the upper state of the clock transition be a long-lived 
metastable state. One type of clock is implemented us- 
ing a group of cold atoms trapped in an optical lattice. 
The second consists of a single laser cooled ion. With ex- 
tremely low systematic perturbations and better stability 
and accuracy, such optical frequency standards should 
exceed the performance of the existing Cs standard. A 
commonly quoted target for the new generation of op- 
tical frequency st anda r ds is a fr actional uncertainty of 

Ai//i/ = icr 18 d, HH, HM . 

There are two main interconnecting areas of theoretic 
atomic clock research: prediction of atomic properties 
required for new clock proposals and determination of 
quantities contributing to the uncertainty budget. New 
clock proposals require estimates of the atomic properties 
for details of the proposals (transition rates, lifetimes, 
branching ratios, magic wavelengths, scattering rates, 



etc.) and evaluation of the systematic shifts (Zeeman 
shift, electric quadruple shift, blackbody radiation shift, 
AC Stark shifts due to various laser fields, etc.). While 
a large fraction of these quantities may be eventually 
measured, lack of knowledge of some of these properties 
may delay experimental realization of new proposals. In 
the case of well-developed proposals, one of the main un- 
certainty issues is the blackbody radiation (BBR) shift. 
The operation of atomic clocks is generally carried out at 
room temperature, whereas the definition of the second 
refers to the clock transition in an atom at absolute zero. 
This implies that the clock transition frequency should 
be corrected for effects of finite temperature, of which the 
leading contributor is the blackbody radiation shift. The 
BBR shift is looming as a major component in the uncer- 
tainty budget of the optical frequency standards. Table 
IXVIII shows the fractional u ncer tainty budget for a 87 Sr 
optical frequency standard [347]. The BBR shift is by 
far the largest source of uncertainty in the uncertainty 
budget. It is noteworthy that the second largest source 
of uncertainty is the AC stark shift caused by the optical 
lattice (this estimate did not take into account possi- 
ble corrections due to Ml and E2 multipoles caused by 
spatial inhomogeneities of the lattice field). Experimen- 
tal measurements of BBR shifts are difficult and high- 
precision theoretical calculations are presently needed. 



1. Black body radiation shifts 

The BBR shift is the AC Stark shift resulting from 
the ambient blackbody radiation field surrounding the 
atom. The BBR energy shift of an atomic state can be 
approximately calculated as [57j 

AE = --l(a 7 r) 3 a (0)T 4 (l + v ) , (66) 
15 

where a is the fine structure constant. The static scalar 
polarizability ao(0) and energy shift AE in Eq. (1551) are 
in atomic units. In this expression, the temperature in 
K is multiplied by 3.166 8153 x 1(T 6 . This is converted 
to Hz by multiplying by 6.579684 x 10 15 . The factor r\ is 
a correction factor that allows for the frequency depen- 
dence of the polari zabil ity when the blackbody integral is 
performed [111 ll3L 12351 ] . The factor 77, referred to as th e 
dynamic shift, is most conveniently written as [l3l |235| 



4Q7r 2 r 2 a/ 



(67) 



The dynamic shift is largest when the excitation energies 
of the states that make the largest contribution to the 
polarizabilities are small. The dynamic shift is largest 
for strontium and increases the BBR shift by 2.7% [13 ] . 

Under most circumstances, the energy shift of an 
atomic level by a radiation field is dominated by the 
dipole component. However, other multipoles might 
make a contribution when the atomic level is part of 
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TABLE XVI: The blackbody radiation shifts for a number of proposed optical frequency standards. The polarizability differ- 
ence, 5a, is negative when the upper state polarizability is smaller than the lower state polarizability. A negative polarizability 
difference means the frequency shift is positive. All BBR shifts are evaluated at 300 K and values that include the dynamic 
shifts are indicated with an asterisk (*). Linewidths are converted from lifetimes, r using <5^ na t = l/(27rr); natural linewidths 
are given for fermionic isotopes for the ns 2 — nsnp clock transitions. Uncertainties in the last digits are given in parentheses. 
References are given in square brackets. The composite CI calculation for Yb + is a hybrid calculation that used CI to explicitly 
allow for core excitations but also included core polarization using a semi-empirical core polarization potential. 



Transition 


u (xlO 15 Hz) 


A 7A. (Hz) 

■ — * •* n a, l V / 


Aq fan) 


A/^rp.r iHzl 


|A^ B R| xl0 15 


Approach 


Ca+(4si /2 -3d 5 / 2 ) 


0.411 [3261 


0.14 [9] 


-44.1(1.5) 


0.38(1) 


0.925 


RLCCSDT [1921 








-42.8 


0.369 


0.895 


CICP [190] 


Sr+(5s 1/2 -4d 5/2 ) 


0.445 [327] 


0.4 [9J 


-29.3(1.1) 


0.250(9)* 


0.562* 


RLCCSDT [1931 


Hg+(5d lu 6s - 5d^ /2 6s 2 ) 


1.06 [3281 


1.8 [9] 








Cryogenic 


Yb+(4/ 14 6s-4/ 13 6s 2 F 7/2 ) 


0.642 [3291 


~ 10" 9 [9J 


11.7 


-0.101 


0.16 


/-sum composite CI [329-331 








6.9(1.8) 


-0.057(14) 


0.089 


/-sum (Lifetimes) [329, 3311 


Yb+(4/ 14 6s-4/ 14 5d 2 D 3/2 ) 


0.688 [332] 


3.1 [9] 


42(8) 


-0.36(7) 


0.53(10) 


Expt. [3321 


Al+(3s 2 1 S - 3s3p 3 PS) 


1.12 [333] 


0.008 [9] 


0.483 


-0.0042(32) 


0.004(3) 


CICP [2351 










-0.008(3) 


0.007(3) 


Expt. [9JJ 


In+(5s 2 1 S-5s5p 3 Pg) 


1.27 [334. 3351 


0.8 [9] 


< 30.7 


> -0.264 


< 0.20 


Theory, using Aa(Cd) 


Mg(3s 2 1 S - 3s3p 3 P °) 


0.655 [336J 


0.00014 [108J 


29.9(7) 


-0.258(7)* 


0.394(11) 


RCI+MBPT [13J 








30.1 


-0.259 


0.395 


CICP [232] 


Mg(3s 2 1 S - 3s3p 3 P?) 


0.656 [337] 


57 [108] 


30.1 


-0.259 


0.394 


CICP [232] 


Ca(4s 2 1 S -4s4p 3 P °) 


0.454 [33J 


0.0005 [108] 


133.2(2.0) 


-1.171(17)* 


2.58(4) 


RCI+MBPT [13] 








135.9 






CICP [234] 


Sr(5s 2 1 S - 5s5p 3 PS) 


0.429 [338J 


0.0014 [108] 


261.1(3.6) 


-2.354(32)* 


5.49(7) 


RCI+MBPT [13] 


Yb(6s 2 1 S-6s6p 3 PS) 


0.518 [339] 


0.008 [340J 


155(15) 


-1.34(13)* 


2.6(3) 


RCI-MBPT [13] 








161(15) 


-1.39(13) 


2.7(3) 


Hybrid RCI+MBPT [143] 


Zn(4s 2 1 S-4s4p 3 PS) 


0.969 [336J 


0.0025 [341] 


29.57 


-0.255 


0.263 


RCICP [113] 


Cd(5s 2 1 S-5s5p 3 PS) 


0.903 [33J 


~ 10 -2 


30.66 


-0.264 


0.292 


RCICP [113] 


Hg(6s 2 1 S-6s6p 3 P °) 


1.13 [145] 


0.11 [3421 


21.0 


-0.181 


0.160 


RCI+MBPT [145] 








24.00 


-0.207 


0.183 


RCICP [113] 



a spin-orbit multiplet [13j. The nsnp 3 Pj levels of the 
alkaline-earth atoms have relatively small energy split- 
tings. The frequency shift due to magnetic dipole (Ml) 
transitions could become important at the 10~ 18 level of 
accuracy. The Ml frequency shift for Sr has been esti- 
mated at 2.4 x 1(T 5 Hz [13}. The frequency shifts for 
other alkaline earths can be estimated using the approxi- 
mate result 5v x « 6v Sl 6E Sl ( 3 P? - 3 P °) /6E X ( 3 P? - 3 P °) 
since the magnetic dipole matrix elements between the 
two members of the triplet show little variation between 
different species. 

Table IXVll lists the frequencies, linewidths, and black- 
body radiation shifts for a number o f potent i al optical 
fre quency standards from Refs. [ p , [l ^ , j^, Jl ^ [Ti a fT43l 
145 
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342J, |348|- The polarizability difference, Aa, is negative 
when the upper state polarizability is smaller than the 
lower state polarizability. All BBR shifts are evaluated 
at 300 K. Linewidths are converted from lifetimes, t us- 
ing <5z/ nat = l/(27rr); natural linewidths are given for 
fermionic isotopes for the ns 2 — nsnp clock transitions. 
It is immediately apparent that the proposed ion clocks 



generally have smaller polarizability differences than the 
electrically neutral atoms in the lattice clocks. 

All of the proposed frequency standards, with the ex- 
ception of the Al+(3s 2 -3s3p 3 P °) transition have T = 300 
K fractional shifts of 10 -16 or higher. The frequency 
shifts of the neutrals are generally larger than the singly 
charged ions. 



2. Optical lattice clocks 

Many of the issues that impact on the optimal choice 
for an optical frequency standard are present in the 
proposed strontium 1 S'- 3 Pq optical frequency standard 
[338l |349| . While Sr might be a desirable atom from the 
perspective of practical experimentation, it is the most 
susceptible to BBR shifts since the polarizability differ- 
ence between the two states is 259 .8 a.u., giving a BBR 
shift of 2.35 Hz at T = 300 K 



141]. 



Assuming that the polarizability difference can be de- 
termined to 0.2% accuracy, the resulting BBR uncer- 
tainty would be 0.0047 Hz which corresponds to a frac- 
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TABLE XVII: Fractional uncertainty budget for the 87 Sr 
atomic frequency standard |347t |. The BBR shifts are evalu- 
ated at 296 K. Corrections that include knowledge of polar- 
izabilities are preceded by an asterisk (*). 



Effect 


Correction 


Uncertainty 




(xlO 16 ) 


(xlO 16 ) 


^Lattice Stark shifts 


-6.5 


0.5 


^Lattice hyperpolarizability 






Stark shifts 


-0.2 


0.2 


*BBR shifts 


52.1 


1.0 


*Probe laser Stark shifts 


0.2 


0.1 


1st order Zeeman 


0.2 


0.2 


2nd order Zeeman 


0.2 


0.02 


Collisional shift 


8.9 


0.8 


Line pulling 


0.0 


0.2 


Servo error shift 


8.9 


0.5 


2nd order Doppler shift 


8.9 


< 0.01 


Totals 


54.9 


1.5 



tional uncertainty of 1 x 10~ 17 . Achieving such a level 
of precision requires experimental determination of the 
five most important transitions in the oscillator strength 
sum-over-states to a precision of 0.1% [141 1. 

Another problem associated with large polarizability 
differences is the enhanced sensitivity with respect to 
variations in temperature. A 1.0 K uncertainty in the 
temperature at 300 K would lead to a frequency uncer- 
tainty of Ai> = 0.064 Hz. The large BBR shift makes a Sr 
standard particularly sensitive to an imprecisely known 
temperature. These problems can be reduced by run- 
ning the clock at lower temperatures. For example, the 
BBR uncertainties stated above can be reduced in size 
by a factor of more than 200 by maintaining the clock at 
liquid nitrogen temperatures. 

Sensitivity to BBR fields has resulted in a proposal 
that mercury would be a sup erior candidate for an op- 
tical frequency standard 145| despite the inconvenience 
of much shorter optical lattice wavelengths. Cadmium 
and zinc have also been ide ntified as candidates with re- 
duced BBR shifts [m|34l|. The drawback of the group 
IIB atoms are the greater uncertainties in the determina- 
tion of the polarizabilities. The underlying (nd) 10 shell 
of the group IIB atoms implies large core polarizabili- 
ties, stronger valence-core correlations and valence ex- 
pectation values that are slower to converge. In addi- 
tion, the resonant oscillator strength for these atoms is 
about 1.4, as opposed to 1.7-1.8 for the group II atoms. 
Consequently the use of a high precision resonant transi- 
tion matrix element from a photo-association experiment 
would do less to minimize the uncertainty than in a group 
II atom. 

Ytterbium has al so be en the subject of increased ex- 
perimental interest [3391 ]. This system also suffers from 
the drawback that it has a large polarizability. Further- 



more, a first-principles calculation of the polarizability 
to a guaranteed accuracy of even 1% is a very difficult 
proposition. The most weakly bound core shell is the 
(4,n 14 shell and the Yb 2+ polarizability is ^9 a.u. 

The atoms that have so far been used in most experi- 
ments, are those that are amenable to cooling and trap- 
ping. The lighter group II atoms, Be and Mg, have the 
disadvantage that they are difficult to c ool, but have the 
advantage of much smaller BBR shifts [3371 ] . Further, it 
would be easier to compensate for the effect of the BBR 
shift in Be and Mg than in most other atoms. Besides 
having smaller shifts, these are relatively light atoms with 
small core polarizabilities, so the uncertainties associated 
with any calculation will be smaller than those of other 
lattice clock. These considerations apply most strongly 
to beryllium. In this case, the polarizability difference of 
the clock atom states is only 1.8 a.u. (Table |XI|) . Beryl- 
lium has only four electrons, so calculations with ECGs 
are possible and should ultimately be able to achieve a 
precision approaching 0.01 a.u. 

The dynamic correction to the BBR shift makes a fi- 
nite contrib ution when the precision reaches the 10~ 18 
Hz level [13, Il93l | but should not lead to a significant in- 
crease in the BBR shift uncertainty. The sum rule for 
evaluation of S(— 4) is more strongly dominated by a few 
major transitions than ao and the relative uncertainty in 
S(— 4) will not be any larger than that of evo- Further, 
the dynamic contribution will be small so the need for a 
precise evaluation is reduced. 

One recent complication has been the recent realiza- 
tion that higher order multipoles could have an impact 
upon the magic wavelength. The inhomogenous spatial 
distributions of the electric and magnetic fields in the 
standing wave patterns that define the lattice can lead 
to energy shifts in the atomic vibrational motion [350j ]. 
This requires the definition of a motion insensitive magic 
wavelength which requires knowledge of the frequency 
dependent el ectri c quadrupole and magnetic dipole po- 
larizabilities [35 ll. 



3. Ion clocks 

Ion state polarizabilities are generally smaller than 
those for neutral atoms because the electrons are more 
tightly bound. None of the ion clocks have polarizability 
differences that exceed 50 a.u.. 

The Al+(3s 2 -3s3p 3 P °) transition has the smallest 
BBR shift of any ion clock due to the fortuitous near 
equality of polarizabilities of the two states in the clock 
transition. The CICP BBR shift is only -0.0042(32) Hz 
[235l | while experiment gave —0.008(3) Hz [97|. How- 
ever, the technical requirements for construction of an 
Al + clock are much more demanding since the cloc k tran- 
sition and cooling laser are in the ultraviolet [333]. 

While the Ca + system is monovalent, calculation of its 
polarizabilities using RLCCSDT method is a more diffi- 
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TABLE XVIII: Summary of the recent theoretical calculations of the Stark shift coefficient k in 10~ 10 Hz/(V/m) 2 and the 
BBR radiation shift parameter ft for transitions between the ground hyperfine states and comparison with experiment. All 
BBR shifts are evaluated at 300 K. Uncertainties in the last digits are given in parentheses. References are given in square 
brackets. 



Atom 


Transition 


ivietnoa 




R 
P 




ZS \r — Z <rr r — 1J 




— U.U0oz4 


n c^m 7vi n — 1^ 
— U.oUl / X 1U 






H/Xpt. loOzj 


— U.Uol(2 ) 




in a 


OS yr — Z <r-r r — L) 


pt r i r i Q"n , T tq^qi 

rLLULjlJ 1 loOoJ 


— U. IzoO 


n pini Qvi n — 14 
— u.ouiy x iu 






H/Xpt. ^jDzJ 


— U.lz4^o j 




39 is 
K 


4s (F — 2 o F — 1) RLCCSDT 1176] 


—0.0746 


—1.118x10 








H/Xpt. loozj 


— U.U < n^j 




Rb 


os (r = 2 <-> _r = 1) 


KLCCbDl 13541 


— 1.272* 


— 1.2o7 X 1U 






RC1+MBP1 11301 


— 1.24(1) 


— 1.26(1) x 10 






iixpt. loozj 


— 1.2,3(0) 




133 Cs 




RLCCSDT [355] 


-2.271(8) 


-1.710(6)xl0 _14 






Theory, PTSCI [129] 


-2.26(2) 


-1.70(2) xlO -14 






Expt. [3561 


-2.271(4) 


-1.710(3)xl0- 14 






Expt. [357] 


-2.05(4) 


-1.54(4) xlO" 14 


137 Ba+ 


6s (F = 2 o F = 1) 


PTSCI [1301 


-0.284(3) 


-0.245(2) xl0~ 14 


171 Yb+ 


6s (F = 1 o F = 0) 


RMBPT3 [1241 


-0.1796 


-0.0983xl0~ 14 






Theory, PTSCI [1301 


-0.171(9) 


-0.094(5) xlO" 14 


199 Hg+ 


6s (F = 1 <-> F = 0) 


PTSCI [130] 


-0.060(3) 


-0.0102(5)xl0~ 14 



^Preliminary value 



cult proposition than for the iso-electronic neutral potas- 
sium [192]. The difficulties lie in the determination of the 
3d state polarizability. First, the 3d state is quite com- 
pact and its charge distribution does perturb the charge 
distribution of the 3s 2 3p 6 core. This leads to a more 
slowly convergent perturbation theory or CI expansion. 
Second, about 30% of the polarizability comes from the 
3d — > nf excitations. The sum-over-states in this case is 
not dominated by a single transition, so discrete excita- 
tions up to the 12/ have to be included. Furthermore, the 
continuum contribution is significant. Including states up 
to 12/ means a much larger B-spline basis needs to be 
used, which in turn makes the calculation more exacting. 
Similar cons iderations impact the BBR shift calculation 
for Sr+ pa) . 

The relative uncertainties associated with the deter- 
mination of the BBR shifts for the two proposed Yb + 
standards are also large due to the underlying 4/ 14 core. 
This is partly mitigated by the small size of the BBR 
shifts. The BBR shift for In + was determined by assum- 
ing the polarizability difference would be smaller than 
that of cadmium. This estimate will be an overestimate 
since the In+ ion will have smaller polarizabilities than 
cadmium. 



4- Experimental possibilities 

So far discussions have focused largely on theory-based 
approaches to the determination of the relevant polariz- 



abilities. However, experimental avenues do exist. For 
example, the polarizability of the Si 2+ ground state, an 
ion iso-electronic with Al + has been determined by RE- 
SIS to an accuracy of better than 0.1%. A RESIS ex- 
periment on Al + should be able to achieve a similar pre- 
cision. Similarly, a RESIS experiment should be able to 
determine the In + ground state polarizability to an accu- 
racy of 0.1%. However, an improved theoretical analysis 
would be needed to get RESIS polarizabilities for Ca + , 
Sr + and Yb + . Application of RESIS to excited parent 
ions also remains a challenge. 



The actual knowledge of the ground and excited state 
polarizabilities is mainly important because it enables 
the determination of the BBR Stark shift. Direct Stark 
shift experiments, on the other hand, might ultimately 
give the most accurate polarizability differences. Table 
IVIII shows that experiments on the ns-npi/2 transitions 
of the alkali atoms have yielded polarizability differences 
with uncertainties less than 0.1 a.u.. 



Photo-association (PA) experiment lifetimes have been 
utilized in estimating polarizabilities with sub-1% preci- 
sion. However, PA spectroscopy has never been applied 
to measure transitions to an excited state, and excited 
state polarizabilities have significant contributions from 
more than one transition. 
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5. BBR shifts in microwave frequency standards 

A BBR shift also exists for the different hyperffne 
states involved in microwave frequency standards. In the 
case of the optical transitions, the lowest (second) order 
polarizabilities of the clock states are different. In the 
case of the ground-state hyperfine microwave frequency 
standards, the lowest (second) order polarizabilities of 
the clock states are identical and the lowest-order BBR 
shift vanishes. To evaluate the BBR shift, third-order 
-F-dependent polarizabilities must be calculated. 

The third-order F-dependent (F is the angular mo- 
mentum of the hyp erfine state) static polarizability, (Xf 
can be written 355] 



a P = Agifx n (2T + C + R), 



(68) 



where A is an angular coefficient, G/ is the nuclear gy- 
romagnetic ratio, and fi n is the nuclear magneton. The 
quantities T, C and R arise from third-order perturba- 
tion theory and typically involve two electric-dipole ma- 
trix elements (i||£)||j) and a one matrix element involving 
the magnetic hype rfine operator For example, term 
T is given by [355j 

T= J2 i HI^IMHPIWHIT«IH 



(E m — E v ) (E n — E v ) 



Here, A\ is the angular coefficient and sums over m, n 
run over all possible states allowed by the selection rules. 

The BBR shift at room temperature effecting the Cs 
microwave frequency standard has b een calculated to 
high accuracy (0.35% and 1%) in Refs. jl29l . l355| . respec- 
tively, implying a 6 x 10 -17 fractional uncertainty. These 
calcu lations are in agreement with a 0.2% measurement 

A summary of recent theoretical calculations [124 

mi Ezl, I353M355I 1 of the Stark shift coefficient k in 

10~ 10 Hz/(V/m) 2 and the BBR radiation shift parame- 
ter (3 for transitions between the gr ound hyp erfine states 
and comparison with experiment [352l |356| is given in 
Table IXVIII1 All BBR shifts are evaluated at 300 K. 
The Stark coefficient k is defined as 



5v = kE 2 , 



(69) 



where 8v is the frequency shift in the static electric field. 
The Stark coefficient for the transition between states F 
and / is related to the polarizability as 



k = ~[a (F) - a (I)}. 



(70) 



The parameter /3 of the relative temperature-dependent 
BBR shift of the microwave frequency standard is defined 
as 



Su 



(71) 



where To is generally taken to be room temperature, 
300-ftT, e parameterizes the lowest-order (in T) contribu- 
tion to the dynamic correction 77 in Eq. (|66j) . and is 
clock transition frequency. The parameter j3 is calcu- 
lated directly from the Stark-shift coefficient k defined 
by Eqs. (pITTO]) as 

k 

P=— (831.9 V/m) 2 . (72) 



D. Long-range interatomic potentials 

The long-range dispersion interaction between two 
spherically symmetric atoms has the form 



C 



TT W U ^ O —'10 

disp ~~&~R*~R^ 



(73) 



where the C n coefficients are called the dispersion coef- 
ficients. The calculation of the dispersion interaction is 
closely related to polarizability calculations. For exam- 
ple, the Cq parameter for two atoms, a and 6, in states 
m and n, can be evaluated using the oscillator strengths 
as 



Co 



E 



fmifi 



"J 



AE mi AE n j(AE mi + AE n j) 



(74) 



The equation is reminiscent of Eq. (|10l) and any calcu- 
lation using Eq. (|74[) automatically generates the nec- 
essary information to generate the dipole polarizability. 
The dispersion coefficients can also be directly evaluated 
from the polarizability at imaginary frequencies as 



Co = - 



a a>0 (tw) a bj0 (iuj) duj. 



(75) 



The polarizability of state n at imaginary frequencies is 
written 



M = E 



fn 



(76) 



Equation (I73[) and subsequent expressions given by 
Eqs. (174")) , or ([751 are the best way to evaluate long-range 
atom-atom interactions. Orthodox quantum chemistry 
techniques are not well suited to determining the very 
small energies of the long range potential. 

The importance of a good description of the long- 
range atom-atom interaction increases at very low ener- 
gies. Determination of the dissociation energy for many 
molecules often involves an extrapolation from the rovi- 
brational energy levels of the highest vibrational states 
|363j |. This has been accomplished i n th e semi-classical 
(WKB) LeRoy-Bernstein procedure [363]. Similarly, the 
determination of the scattering length in cold-atom col- 
lisions often requires knowledge of the dispersion param- 
eters [364] , 
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TABLE XIX: The lowest order dispersion coefficient, Cq for homo-nuclear atom-atom pairs. The Hybrid-RLCCSD replaces 
the calculated matrix element for the resonance transition with an experimental value. 





00 Li 


Na 




Rb 


V b 


Fr 


Hylleraas [118] 


1393.42(5) 












Model Potential [358J 


1388 


1472 


3813 


4426 


6331 




CICP [105] 


1394.6 


1561 


3905 


4635 






RLCCSD [56] 




1564 


3867 


4628 


6899 


5174 


Hybrid-RLCCSD [56. 3591 


1390(2) 


1556(4) 


3897(15) 


4691(23) 


6851(74) 


5256(89) 


Expt. 






3921 [3601 


4698(4) [3611 


6877(24) [74] 




Expt. 










6860(25) [362] 





Better information about the specific values of the dis- 
persion coefficients for many atoms has become available 
primarily because of the importance of such data for the 
field of cold-atom physics. There have been the near ex- 
act non-relativistic calculation by Yan an d co-workers on 
H, H e and Li using Hylleraas basis sets jll7l Ill8l . 13651 — 
13691 ] . An important series of calculations on the ground 
and excited states of the alkal i atoms w ere reported by 
Marinescu and co-workers |358l . l370l - l373| . However, these 
calculations were performed with a model potential ap- 
proach that omitted some dynamical features (e.g. tran- 
sitions from the core) that should be included. Later 
calculations with sem i -empirica l Hamiltonians by Mitroy 

and co-workers OH [ml IM HHH iil and rl- 

CCSD/RCI+MBPT calculatio ns by Derevianko and co- 
workers [Hll, HE H3, Hz! Effif should be preferred 
since the underlying atomic structure descriptions are su- 
perior. These calculations encompass both the alkali and 
alkaline-earth atoms. Table IXIXI shows that CICP and 
RLCCSDT calculations of Cq for homo-nuclear pairs of 
alkali atoms agree at the 1% level. This ag reem ent ex- 
tends to hetero-nuclear pairs of alkali atoms [105| and to 
alkaline-earth atoms 



ity experiment admits other interpretations. Taking the 
polarizability and diamagnetic susceptibility as known 
quantities, the refractive index experiment yields a value 
for the universal gas constant, R = 8.314487(76), which 
is not far removed in prec ision from the recommended 
value of 8.314472(15) [373]. Boltzmann's constant, the 
definition of the mol and the universal gas constant are 
all inter-related through the identity, R = IcbNa- 



F. Atomic transition rate determinations 

The sum-over-states approach described in Sec- 
tion [TTX2] is generally used to determine the polarizabil- 
ities from calculated or experimental oscillator strengths 
or El matrix elements. It is possible to reverse the pro- 
cess for systems which have a precisely known polariz- 
ability that is dominated by a single strong tran sition. 
A good example occurs for the cesium atom [3801 ] where 



the dipole polarizability 74] and line strength ratio [381 ] 
have been measured to high accuracy. 

The ground state static polarizability ao can be writ- 
ten as 



E. Thermometry and other macroscopic standards 

The present definition of temperature is based on the 
triple point of water which is set to 273.16 K. An alter- 
native approach would be to fix the Boltzmann constant, 
ks and then measure the thermometric properties of a 
substance which depend on the product ksT. At present, 
the best estimate of the Boltzmann constant was deter- 
mined by the speed of sound in helium gas. Acoustic 
gas thermometry (AGT) has resulted in a value of ks 
accurate to 1.8 ppm [3771 1371 ] . 

The speed of sound is not the only thermometric prop- 
erty that can be used to determine ks- Two other proper- 
ties are the dielectric constant fo r helium gas and the re- 
fractive index for helium gas [H, 13781 . l379j . The most re- 
cent refractive index experiment using a microwave cav- 
ity [63| has given the dipole polarizability to an accu- 
racy of 9.3 ppm. If the 4 He polarizability is taken as 
a known quantity from theory, then the microwave cav- 



(77) 



where ae p is the contribution of the resonance excitations 
to the polarizability, i.e. from 6s — 6pi/2 and 6s — 6p 3 / 2 
transitions, and a' v includes contributions from all other 
excited states . Rearranging and expressing ag p in terms 
of the 6s — > 6pi /2 line strength gives 



Sqs— 6pi/2 



"0 - Ot v ~ "c 



1 



R 



(78) 



3A£, 



6^-6pi/ 2 



3AE 



6s-6p 3/2 



The factor R is the ratio of the line strengths of the spin- 
orbit doublet. Using a = 401.0(6) 74], R = 1.9809(9) 
[38lj . ionic core polarizability of 15.644(4) jTT^ that 
needs to be corrected for the presence of the valence 
electron by the term a cv — —0.72, a' v — 1.81 yields 



6s-6pi/ 2 



20.308(42) and 5, 



'6a-6 P3/2 = 40.227(84). The 
corresponding values for the reduced matrix elements in 
atomic units are 4.510(4) and 6.347(5) for the 6s — 6pi/2 
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and 6s — 6j<3/2 transitions respectively. The uncertainties 
of these values are dominated by the uncertainty in the 
experimental value of o.q. 

A similar approach has been used to determine the 
multiplet strengths for the resonance transitions in Mg + , 
Si 3+ J81 and Si 2+ 88] from RESIS experimental data 
[198L l239l ]. 

Stark shifts for the ns-npi/ 2 transition [2 Oil ] have also 
been used to derive estimates for the S(npi/ 2 — (^— 1)^3/2) 
line strengths with a precision of about 1% for potassium 
and rubidium [301 ] . This analysis relied on the result that 
80-90% of the npi/ 2 polarizability comes from the exci- 
tation to the (n — l)d 3 /2 state. These values were also 
used to determine the magic wavelengths for the np — ns 
transitions in these alkali atoms (29| . Such determination 
of matrix elements permitted benchmark comparisons of 
theory and experiment [30j . 

The procedures described above also permit the cross- 
checking of results from completely different types of ex- 
periment. The domination of the 6p Cs scalar polar- 
izabilities by the 5d — 6p dipole matrix elements facil- 
itated an exacting consiste ncy check of the 5d lifetime 
with 6p polarizability data 169]. In that work, 5d — 6p 
matrix elements obtained from experimental Stark shift 
data were compared with the values extracted from the 
5d lifetimes. The experimental measurements of the 5d 
lifetime and 6p scalar polarizabilities were found to be 
inconsistent within the uncertainties quoted by the ex- 
periment al gro ups 169]. Theoretical RLCCSDT matrix 
elements [169j were found to be in agreement with the 
Stark shift experiments but not with the lifetime mea- 
surements. 



VIII. CONCLUSIONS 

The advent of cold-atom physics owes its existence to 
the ability to manipulate groups of atoms with electro- 
magnetic fields. Consequently, many topics in the area 
of field-atom interactions have recently been the subject 
of considerable interest and heightened importance. This 
applies to a quantity like the dipole polarizability which 
governs the first-order response of an atom to an applied 
electric field and the preceding few years have seen many 
calculations of atomic polarizabilities for a variety of sys- 
tems. 

The aim of the present review has been to provide 
a reasonably comprehensive treatment of polarizability 
related issues as they relate to topics of contemporary 
importance. However, our treatment is not exhaustive. 
The polarizabilities of many atoms such as the halogens 
have been omitted. The reader is referred to the broader 
treatment in [2l[. Similarly, the treatment of DC and 
AC Stark shift data is better described as selective as 
opposed to exhaustive. 

Part of the motivation for this review has been the im- 
portance in developing new atomic based standards of 
time 344] . and corresponding need for precise knowledge 



of the blackbody radiation shifts. The primary require- 
ment for the BBR application is for polarizabilities and 
Stark shifts to be known with a precision of 0.1% or bet- 
ter. Much of the existing body of experimental data is an 
order of magnitude less precise. Direct measurements of 
clock transition Stark shifts would be helpful in reducing 
the BBR shift uncertainties. 

One area where theory might be useful in this endeavor 
would be in the development of atom based polarizability 
standards. Such a standard is already in existence for he- 
lium where theoretical and experimental polarizabilities 
have uncertainties of 0.17 ppm and 9.1 ppm [63|, respec- 
tively. These results are not relevant to atomic clock 
research and another atom needs to serve as a standard. 
Hylleraas calculations on lithium could yet serve to pro- 
vide a theoretical reference point for Stark shift experi- 
ments. At the moment the uncertainties in t he best cal- 
culation and best experiment are 0.11% [l T0| | and 0.07% 
9l|. A better treatment of rclativistic effects should re- 
sult in uncertainty in the Hylleraas calculation decreasing 
to the 0.01% level of precision. 

One possible avenue for improvement could be in the 
development of hybrid theoretical approaches combin- 
ing the best features of different methods. For example, 
orbital-based approaches cannot match the extreme ac- 
curacies achievable with correlated basis sets. Direct in- 
corporation of the the Dirac Hamiltonian in orbital-based 
calculations is now relatively routine, but this is not the 
case for calculations with correlated basis sets. Perhaps, 
comparisons of correlated basis calculations with non- 
relativistic orbital-based calculations and with relativis- 
tic orbital-based calculations could be used to estimate 
relativistic corrections to Hylleraas calculations or corre- 
lation corrections to orbital-based calculations. 

It is likely that the determination of polarizabilities will 
become increasingly important in the future. As experi- 
ments become capable of greater precision, it will become 
necessary to make more detailed corrections of the effects 
of electromagnetic fields that are used for manipulation 
and investigation of atoms. 
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